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Abstract 

Consider a system of N bosons in three dimensions interacting via a repulsive short range 
pair potential N 2 V(N(xi — Xj)), where x = {x\, . . . ,xm) denotes the positions of the particles. 
Let Hpf denote the Hamiltonian of the system and let if>N,t be the solution to the Schrodinger 
equation. Suppose that the initial data ipN,o satisfies the energy condition 

(iPn,o,H^n,o) <C k N k 

for k — 1,2,.... We also assume that the fc-particle density matrices of the initial state are 
asymptotically factorized as N — > oo. We prove that the fc-particle density matrices of i/jN,t 
are also asymptotically factorized and the one particle orbital wave function solves the Gross- 
Pitaevskii equation, a cubic non-linear Schrodinger equation with the coupling constant given by 
the scattering length of the potential V. We also prove the same conclusion if the energy condition 
holds only for k — 1 but the factorization of ^>jv,o is assumed in a stronger sense. 

1 Introduction 



Bose-Einstein condensation states that at a very low temperature Bose systems with a pair interaction 
exhibit a collective mode, the Bose-Einstein condensate. If one neglects the interaction and treats 
all bosons as independent particles, Bose-Einstein condensation is a simple exercise |15j . The many- 
body effects were traditionally treated by the Bogoliubov approximation, which postulates that the 
ratio between the non-condensate and the condensate is small. The coupling constant a/8ir obtained 
by the Bogoliubov approximation is the semiclassical approximation of the the scattering length 
ao of the pair potential. To recover the scattering length, one needs to perform a higher order 
diagrammatic re-summation, a procedure that yet lacks mathematical rigor for interacting systems. 

Gross |12t [T3] and Pitaevskii [20] proposed to model the many-body effects by a nonlinear on-site 
self interaction of a complex order parameter (the "condensate wave function" ) . The strength of the 
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nonlinear interaction in this model is given by the scattering length ao- The Gross-Pitaevskii (GP) 
equation is given by 
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where £ is the Gross-Pitaevskii energy functional and a = 8iraQ. The Gross-Pitaevskii equation 
is a phenomenological mean field type equation and its validity needs to be established from the 
Schrodinger equation with the Hamiltonian given by the pair interaction. 

The first rigorous result concerning the many-body effects of the Bose gas was Dyson's estimate 
of the ground state energy. Dyson [5] proved the correct leading upper bound to the energy and a 
lower bound off by a factor around 10. Dyson's upper bound was obtained by using trial functions 
with short range two-body correlations. This short scale structure is crucial for the emergence of the 
scattering length and thus for the correct energy. The matching lower bound to the leading order in 
the low density regime was obtained by Lieb and Yngvason [19] . Lieb and Seiringer [IB] later proved 
that the minimizer of the Gross-Pitaevskii energy functional correctly describes the ground state of 
an iV-boson system in the limit N — > oo provided that the length scale of the pair potential is of 
order 1/N. For a review on related results, see |17j . 

The experiments on the Bose-Einstein condensation were conducted by observing the dynamics 
of the condensate when the confining traps are removed. Since the ground state of the system with 
traps will no longer be the ground state without traps, the validity of the Gross-Pitaevskii equation 
for predicting the experimental outcomes asserts that the approximation of the many-body effects 
by a nonlinear on-site self interaction of the order parameter applies to a certain class of excited 
states and their subsequent time evolution as well. 

In this paper, we shall prove that the Gross-Pitaevskii equation actually describes the dynamics of 
a large class of initial states. The allowed initial states include wave functions with the characteristic 
short scale two-body correlation structure of the ground state and also wave functions of product 
form. Notice that product wave functions do not have this characteristic short scale structure, 
nevertheless the GP evolution equation applies to them. It should be noted that our theorems 
concern only the evolution of the one particle density matrix but not its energy. In fact, for product 
initial states, the GP theory is correct on the level of density matrix, but not on the level of the 
energy. We shall discuss this surprising fact in more details in Section [3l 

2 The Main Results 



Recall that the Gross-Pitaevskii energy functional correctly describes the energy in the large N limit 
provided that the scattering length is of order 1/N |18j . We thus choose the interaction potential to 
be 

V N (x) := N 2 V (Nx) = — N 3 V (Nx) . 

This potential can also be viewed as an approximate delta function on scale 1/N with a prefactor 
1 /N which we will interpret as the mean field average. The Hamiltonian of the Bose system is given 
by 



v 



N 



H N :=-J2 A j + J2 V ^ 

j=i j<k 



x k ) , V N (x) := iV 2 V (Nx) . 



(2.1) 



The support of the initial state will not be scaled with N. Thus the density of the system is N 
and the typical inter-particle distance is iV -1 / 3 , which is much bigger than the length scale of the 
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potential. The system is really a dilute gas scaled in such a way that the size of the total system is 
independent of N. 

The dynamics of the system is governed by the Schrodinger equation 

idtipN,t = H N ip N>t (2.2) 

for the wave function tpN jt G L 2 (R 3N ), the subspace of L 2 (R 3N ) consisting of all functions symmetric 
with respect to any permutation of the N particles. We choose ipN,t to have L 2 -norm equal to one, 

UnA = i. 

Instead of describing the system through the wave function, we can describe it by a density 
matrix 7^ G £ 1 (Lg(M 3iV )), where £ 1 (L^(IR 3Ar )) denotes the space of trace class operators on the 
Hilbert space L^(IR 3Ar ). A density matrix is a non-negative trace class operator with trace equal to 
one. For the pure state described by the wave function ip^, the density matrix 7^ = \^n){^n\ is 
the orthogonal projection onto ip^. The time evolution of a density matrix 7jv is determined by the 
Heisenberg equation 

idtlN,t = [HN,lN,t], (2.3) 
where [A, B] = AB — BA is the commutator. 
Introduce the shorthand notation 



x := (x 1 ,x 2 , ■ ■ ■ ,x N ), x fe := (xi, . . . ,x k ), x N _ k := (x k+1 , . . . ,x N ) 

and similarly for the primed variables, x' fc := (x[, . . . , x' k ). For k = 1,...,N, the k-particle re- 
duced density matrix (or k-particle marginal) associated with 7jv,t is the non-negative operator in 
>C 1 (L^(M 3fe )) defined by taking the partial trace of 7^,4 over N — k variables. In other words, the 

(k) 

kernel of *y N t is given by 

7^(x fc ;x' fc ) := y"dx JV _fe7iv,t (xfe^jv-fejx^xjv-fe) . (2.4) 

(k) 

Our normalization implies that Tr j K N ' t = 1 for all k = 1, . . . , N and for every t£l. 

We now define a topology on the density matrices. We denote by C\ = £ 1 (L 2 (M 3fc )) the space of 
trace class operators acting on the Hilbert space L 2 (R 3k ). Moreover, K k = K,(L 2 (R 3k )) will denote 
the space of compact operators acting on L 2 (M. 3k ) equipped with the operator norm, || • \\fc k := || • ||. 
Since C\ = /C£, we can define the weak* topology on C l (L 2 (M. 3k )) , i.e., u n —>■ 00 if and only if for 
every compact operator J on L 2 (R 3k ) we have 

lim Tr J oo n = Tr Jco . (2.5) 

n— >oo 



Throughout the paper we will assume that the unsealed interaction potential, V(x), is a nonneg- 
ative, smooth, spherically symmetric function with a compact support in the ball of radius R, 

supp V C {x e M. 3 : \x\<R}. (2.6) 

With the notation r = |x|, we will sometimes write V(r) for V(x). We define the following dimen- 
sionless quantity to measure the strength of V 

poo 

p := sup r 2 V{r) + / drrV{r) . (2.7) 

r>0 Jo 
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Let / be the zero energy scattering solution associated with V with normalization lim| : j.|_ >00 fix) = 
1. We will write f(x) = 1 — wo(x). By definition, this function satisfies the equation 

A + ^(x)](l-wo(z))=0, (2.8) 
and limui^oo wq(x) = 0. The scattering length ao of V is defined by 

ao := lim u>o(£)|x| . (2-9) 

Since V has a compact support (|2.6p . we have 

/( x ) = l_^ |x|>i2. (2.10) 

From the zero energy equation, we also have the identity 

J dx V(x)(l - w (x)) = 8vra . (2.11) 

By scaling, the scattering length of the potential V/v( x ) is a := ao/N and the zero energy scattering 
equation for the potential Vn is given by 

-A + ^Vn{x) \ (1 — w(x)) = (2.12) 

where w(x) := wq(Nx). Note that w(x) = a/\x\, for \x\ > R/N. 
We can now state our main theorems. 

Theorem 2.1. Suppose V > is a smooth, compactly supported, spherically symmetric potential 
with scattering length oq and assume that p ( defined in (|2.7[) ) is small enough. We consider a family 
of systems described by initial wave functions ipN ^ L 2 S (M? N ) such that 

(iP n ,HUn) <C k N k (2.13) 

for all k > 1. We assume that the marginal densities associated with ipN factorize in the limit 
N — ► oo, i.e. there is a function ip £ L 2 (M 3 ) such that for every k > 1, 



(2-14) 

as N — > oo with respect to the weak* topology of C 1 (L 2 (M 3k )). Then ip E ^(M 3 ), and for every fixed 
k > 1 and t € R, we /iaue 

(2-15) 

with respect to the same topology. Here ipt £ -fT 1 (]R 3 ) is the solution of the nonlinear Gross-Pitaevskii 
equation 

id t <pt = -A<pt + 8ira \<pt\ 2 <Pt (2.16) 

with initial condition pt=o = f- 

Using an approximation argument, we can relax the energy condition (|2,13p . and only assume 
that (ipN, H^iPn) < CN. However, in order to apply our approximation argument, we need to 
assume stronger asymptotic factorization properties on ifj^. 
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Theorem 2.2. Suppose V > is a smooth, compactly supported, spherically symmetric potential 
with scattering length oq and assume that p ( defined in \2. 7| ) ) is small enough. We consider a family 
of systems described by initial wave functions ip^ G L 2 (M 3Ar ) such that 

(iP N ,H Nl p N ) <CN. (2.17) 

We assume asymptotic factorization of ipw in the sense that there exists <p G L 2 (R 3 ) and, for every 
N, and every 1 < k < N, there exists a ^~ k) G L 2 (R 3 ( 7V ~ fc )) with \\^ ~ k) \\ = 1 such that 

UN-V® k ®Z { N~ k) \\^0 (2-18) 
as N — > oo. This implies, in particular that, for every k > 1, 

-rP^\<f)(<p\® k (2.i9) 

as N — > oo with respect to the weak* topology o/£ 1 (L 2 (R 3fc )). Then tp G //-'-(R 3 ), and for every fixed 
k > 1 anc? t£l we /iaue 

7 g - |^}<^r (2-20) 

with respect to the same topology. Here tpt G ^(R 3 ) is i/ie solution of the nonlinear Gross-Pitaevskii 
equation 

id t <pt = -A<Pt + 87ra |^| 2 ^ (2-21) 

with <pt=o = p. 

Both theorems have analogous versions for initial data describing mixed states (that is jn is not 
an orthogonal projection). For example, suppose that 7^ is a family of density matrices satisfying 

Tr H k NlN < C k N k and 7 }f ) -► u® k (2.22) 
where ujq is a one-particle density matrix and 

k 

uf k (x k ;x' k ) = Yl^oix^x'j). 

3=1 

Then for every tGl and k > 1 we have 

7$ - < (2-23) 
where is the solution of the nonlinear Hartree equation 

id t uj t = [-A + 8ira g t ,uJt\ Qt(x) = u t (x; x), uj t=0 = iv (2.24) 
The last equation is equivalent to (|2.16|) if Ut = \ipt}((pt |- 

Lieb and Seiringer |16] have proved that, for pure states, the assumption 

7at -"MM as N 00 

implies automatically (|2.14p for every > 1 (see the argument after Theorem 1 in that paper |3- For 
mixed initial states we still need the second condition in (|2.22p for all k > 1 in order to prove (|2.23|) . 



1 We thank Robert Seiringer for pointing out this result to us. 
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Now we comment on the assumption of asymptotic factorization (|2.18p for the initial data ipN- 
The most natural example that satisfies this condition is the factorized wave function tpjsr(x) = 
1X7=1 f{ x j)- ^ additionally, 99 £ H 1 ^ 3 ), then (|2.17p is also satisfied by the Schwarz and Sobolev 
inequalities. The evolution of ipN is therefore governed by the GP equation according to Theorem 
12.21 This is, however, somewhat surprising because the emergence of the scattering length in the 
GP equation indicates that the wave function has a characteristic short scale correlation structure, 
which is clearly absent in the factorized initial data. We shall discuss this issue in more details in 
Section 

From the physical point of view, however, the product initial wave function is not the most 
relevant one. In real physical experiments, the initial state is prepared by cooling down a trapped 
Bose gas at extremely low temperatures. This state can be modelled by the ground state V , ^ ap °f 
the Hamiltonian 

N N 
j=l i<j 

with a trapping potential V ex t(x) — ► 00 as |x| — > 00. In Appendix [Cl we prove that assumptions 
(12.17j) and (|2.18j) are satisfied for V^ aP • I n other words, Theorem 12.21 can be used to describe the 
evolution of the ground state of -ff^ ap , after the traps are removed (see Corollary lC.ljl . This provides 
a mathematically rigorous analysis of recent experiments in condensed matter physics, where the 
evolution of initially trapped Bose-Einstein condensates is observed. 

In Appendix |Bl we show that Theorem 12.21 can also be applied to a general class of initial data, 
which are in some sense close to the ground state of the Hamiltonian ii^ ap . The ground state of a 
dilute Bose system with interaction potential is believed to be very close to the form 

WNW'-HfiNixi-Xj)), (2.25) 

i<j 

where / = 1 — wq is the zero-energy solution (12.8f) . We remark that Dyson [5] used a different function 
which was not symmetric, but the short distance behavior was the same as in Wjy. An example of 
a family of initial wave functions which have local structure given by W is given by wave functions 
of the type 

N 

il> N (x) = W N (x)Y[<p(x j ) (2.26) 

3=1 

where (p G H 1 ^ 3 ). Due to the factor Wn, this function carries the characteristic short scale 
structure of the ground state. We will prove in Lemma IB.ll that wave functions of the form (|2.26p 
(with correlations cutoff at length scales £ 3> A -1 ) satisfies the assumptions (I2.17P and 112. 18|) . 

Part of Theorem 12.21 was proved in [8] for systems with the pair interaction cut off whenever 

three or more particles are much closer to each other than the mean particle distance, A -1 / 3 . For 

(k) 

this model, it was proved that any limiting point of 7^ satisfies the infinite BBGKY hierarchy 
(see Section [3]) with coupling constant 87rao. The uniqueness of the solution to the hierarchy was 
established in [9]. In the current paper we remove this cutoff and establish the apriori bounds needed 
for the uniqueness theorem in [9]. 

The Hamiltonian (12.11) is a special case of the Hamiltonian 

N N 

Hp,N := - Y, A i + n E NWV ^^ ~ x i)) ( 2 - 27 ) 

j=l i<j 
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introduced in [6] and [9]. In [9] we have proved a version of Theorem 12.21 for < f3 < 1/2 provided 
the initial data is given by a product state V'at(x) = Yij=i f( x j) f° r some cp G if 1 (M 3 ). In this case 
the limiting macroscopic equation was given by 

idttpt = -A(p t + bo\tp t \ 2 tpt , 

with b = J dxV(x). Note that N 3/3 V(N^x) is an approximate delta function on a scale much 
bigger than 0(1/N), the scattering length of j^Vn. This explains why the strength of the on-site 
potential is given by the semiclassical approximation bo of the 87rao- With the techniques used in 
this paper, it is straight-forward to extend the result of [9] to all (3 < 1 with the same coefficient 
&o in the limiting one-body equation provided that p (from (|2.7p ) is small enough. Combining this 
comment with Theorem 12.11 and 12.21 we have shown that the one particle density matrix for the 
iV-body Schrodinger equation with Hamiltonian given by (|2.27p converges to the Gross-Pitaevskii 
equation with coupling constant given by 

if < < 1 

(2.28) 

if p = 1 . 

The case P = is the mean-field case and the limiting one-body equation is the Hartree equation: 

id m = -A<pt + (V * \<p t \ 2 )<Pt ■ (2.29) 

This was established by Hepp |14| for smooth potential. Ginibre and Velo [TTJ considered singular 
potentials but with a specific initial data based on second quantized formalism. Spohn [22] introduced 
a new approach to this problem using the BBGKY hierarchy. Recent progresses on mean-field limit 
of quantum dynamics have been based on the BBGKY hierarchy and we mention only a few: the 
Coulomb potential case OHO], the pseudo-relativistic Hamiltonian with Newtonian interaction [7], 
and the delta function interaction in one dimension by Adami, Bardos, Golse and Teta [1] [2]. In 
next section, we review the BBGKY hierarchy and the two-scale nature of the eigenfunctions of 
interacting Bose systems. 

3 The BBGKY Hierarchy 

(k) 

The time evolution of the density matrices Tjy-i, for k = 1,...,N, is given by a hierarchy of N 
equations, commonly known as the BBGKY hierarchy: 

k k 



7=1 i<j 

k (3-30) 



+ (N - k) Tr fc+ i Vn(xj - x k+ i),T N: 

3=1 



(fe+1) 
t 



for k = 1, . . . , N (we use the convention that 7^ = if k > N). Here Tr^ + i denotes the partial 
trace over the (k + l)-th particle. In particular, the density matrix 7^ t (xi;x' 1 ) satisfies the equation 



id tTN,t(. x r,x' 1 ) = (-A X1 + A x ./)7^(xi;x' 1 ) 

+ (N-1) / dx 2 (V N (x 1 - x 2 ) - V N (x[ - x 2 )) ^ > t (xi,x 2 ;x' 1 ,x 2 ), 



(2) , , (3-31) 
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To close this equation, one needs to assume some relation between jffi t and The simplest 

assumption would be the factorization property, i.e., 

7/?,t0ci, x ^ x '\i 4) = tS(^i; ^1)7^(^2; 4) • ( 3 - 32 ) 

This does not hold for finite N, but it may hold for a limit point 7^ ' of 7)^ as N — > 00, i.e., 

7 i (2) (a;i,a;2;a;i,a; / 2 ) = 7 t (1) (xi; x'^7^ (x 2 ; x' 2 ) . (3.33) 
Under this assumption, 7^ satisfies the limiting equation 

t 7t (1) (^i;4) = ("A,, + A 4 )7 t (1) (x 1 ;xi) + (Q t (n) - Q t (a4)) 7 t (1) (^i;4) (3-34) 

where 

.</> A \-< - '/ /'/>.(; 

If pt{x) is continuous, then Q t is given by 



Q t (x) := lira iV / dyV N (x - y)p t (y), pt(x) = % '(x;x). (3.35) 

iv— >oo 



Qt(x) = b p t (x). 

Thus (|3.34p gives the GP equation with a coupling constant cr = &o instead of <r = 87rao. This 
explains the case if (3 < 1. For /? = 1, we note that &o/8tt is the first Born approximation to the 
scattering length ao and the following inequality holds: 

a <^ = ]- [ \ V(x) dx . (3.36) 

Recall that the ground state of a dilute Bose system with interaction potential Vjv is believed to be 
very close to W(x) (see (|2.25p ). We assume, for the moment, that the ansatz, ^(x) = W(x)0t(x) 
with fa a product function, holds for all time. The reduced density matrices for ifit( x ) satisfy 

l¥\xi,X2;x' 1 ,x 2 ) ~/(jV(xi - x 2 ))f(N(x' 1 - x' 2 ))it 1 \x 1 ] x^jj^fa] x' 2 ) . (3.37) 
Together with (12. lip and the assumption that pt is smooth on scale 1/N, we have 

lim N / dx 2 V N (xi - x 2 )^\{x 1 ,x 2 ;x' l ,x 2 ) = 87ra 7 t (1) (xi; x'^p^xi) . (3.38) 

N^oo J ' 

This formula is valid for \xi — x[\ S> 1/N. We have used that ]im.\ x \^ ao f (x) = 1. For pure states, 
this gives the GP equation with the correct dependence on the scattering length. 

Notice that the correlation in 7^) occurs at the scale 1/N, which vanishes in a weak limit 
and the product relation (|3.33p will hold. However, this short distance correlation shows up in 
the GP equation due to the singular potential NVn(xi — x 2 ). This phenomena occurs for the 
ground state as proved in [18]. Our task is to characterize wave functions with this short scale 
structure and establish it for the time evolved states. The key observation is the following Proposition. 
Recall the assumptions on V from Section [2] and that 1 — w(x) denotes the zero energy solution to 
—A + |Viv (|2.12p . We will use the short notation Wij := w(xi — Xj), Vwij = (Vw)(xi — Xj) (note 
that Viv^ = —Vwji). 
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Proposition 3.1 (ii^-energy estimate). Suppose that p (defined in \2.7\) ) is small enough. Then, 
there exists a universal constant c > such that, for every tp 6 Lg(R^), and for every fixed indices 
i J> h3 = 1) • • • > N, we have 

(V, H 2 N ijj) > (1 - cp)N(N - 1) J (1 - Wij f \ViVj <^| 2 (3.39) 

where 4>ij defined by ip = (1 — u>ij)4>ij. 

If 4>ij is singular when Xi approaches Xj, then VjVj^j cannot be L 2 -integrable. This Proposition 
thus shows that the short distance behavior of any function ip with (ip,H%ifi) < CN 2 is given by 
(1 — w(xi — Xj)) when is near Xj. 

We emphasized the importance of the local structure (1 — w(xi — Xj)) for obtaining the scattering 
length do- While Theorem 12.21 concerns only the one particle density matrix in the weak limit and 
no statement on the local structure is made at all, the validity of the GP equation does suggest the 
existence of this structure. For the initial data (12.260 beginning with this local structure, it simply 
means its preservation by the dynamics. This is indeed the case if the local structure of the initial 
data tp is precise enough so that (ip,H^ip) < CN 2 , see Proposition EH 

For the product initial state, there is no such structure to begin with. Theorem 12.21 thus indicates 
that on some short length scale a local structure similar to (1 — w(x{ — Xj)) forms in a very short time 
which approaches zero in the limit N —* oo. Heuristically, notice that the two particle dynamics is 
described by the operator 

id t - A X1 - A X2 - V N ( Xl - x 2 ) = N 2 [i&r - A Xl - A Xa - V{X X - X 2 )) 

where Xi = Nxi and T = N 2 t are the microscopic coordinates. The small positive time behavior 
of the original wave function on the short length scale is the same as the long time behavior in 
the microscopic coordinates. Clearly, we expect the long time dynamics to be characterized by the 
relaxation to the zero energy solution. This picture, however, is far from rigorous as the true iV-body 
dynamics develops higher order correlations as well. 

On the other hand, the local structure (1 — w(xi — Xj)) cannot be the only singular piece of the 
wave function in positive time for product initial states. A simple calculation shows that the energy 
per particle of a product initial state ^jv(x) = Y[f=i f( x j) ^ s given by 

lim N- 1 (*p N ,H N ip N )= [ dx\Vip(x)\ 2 + ^- [ dx\ip(xt (3.40) 

N^OO J R 3 2 J R 3 

where bo = J V. This is different from the GP energy functional (II. 1|) due to the coupling constant. 
Since the energy is a constant of the motion, this implies that the GP theory does not predict the 
evolution of the energy. If we grant that the local structure (1 — w(xi — Xj)) does form for positive 
time t > 0, the discrepancy in energy suggests that there is some energy on intermediate length scales 
of order N~ a , < a < 1 which is not captured by the GP theory. This excess energy apparently 
does not participate in the evolution of the density matrix on length scale of order one which is the 
only scale that is visible by our weak limit. We do not know if such a picture can be established 
rigorously. 

Notation. We will denote an arbitrary constant by C. In general C can depend on the choice 
of the unsealed potential V. Universal constants, independent of V, will be denoted by c. We write 
f(N) = o(N a ) if there is 5 > such that iV" Q+5 /(iV) -> as N -> oo (unless stated otherwise, 
this convergence does not need to be uniform in the other relevant parameters). We also write 
f(N) <C g(N) if f(N)/g(N) = o(l). Integrations without specified domains are always understood 
on the whole space (M 3 , M. 3k or M. 3N according to the integrand) with the Lebesgue measure. 
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4 Proof of Theorem 12.11 and Theorem 12.2 



In this section we present the main steps of the proofs and we reduce the argument to a sequence of 
key theorems and propositions. These will be proven in the rest of the paper. 

We start with defining the space of density matrices that depend continuously on the time 
parameter with respect to the weak* topology. To use Arzela-Ascoli compactness argument, we 
will need to establish the concept of uniform continuity in this space, thus we have to metrize the 
weak* topology. 

Since Kk is separable, we can fix a dense countable subset of the unit ball of ICk- we denote it 
by {Ji^}i>i £ ICk, with || J^HjCj. < 1 f° r a U i > 1- Using the operators we define the following 
metric on L\: for 7W ,7^ £ C\ we set 

00 



i=l 

Then the topology induced by the metric r\ k and the weak* topology are equivalent on the unit ball 
of C\ (see [2T], Theorem 3.16) and hence on any ball of finite radius as well. In other words, a 

uniformly bounded sequence 7^ £ C\ converges to 7^ £ C\ with respect to the weak* topology, if 

and only if 77^(7]^, 7^)) — > as N — > 00. 

For a fixed T > 0, let C([0,T],££) be the space of functions of t £ [0, T] with values in C\ which 
are continuous with respect to the metric On C([0,T],£^) we define the metric 

%(7 (fe) (-),7 (fc) (-)):= sup % ( 7 ( fc )(t),7( fc )(i)). (4.2) 
te[o,T] 

Finally, we denote by r pro( i the topology on the space © fe>1 ^([0, T], C\) given by the product of the 
topologies generated by the metrics r}k on C([0, T], C\). 



Proof of Theorem \2.1\ The proof is divided in several steps 



Step 1. Compactness ofTN,t = {7w*}fc>i- We set T > and work on the interval i £ [0, T]. 



Negative times can be handled analogously. We will prove in Theorem 16.11 that the sequence T N t = 

{7ivt}fc>i e ©fc>i ^([0j^]>^jfc) i s compact with respect to the product topology r pro d defined above 

(we use the convention that 7^ = if k > N). Moreover, we also prove in Theorem 16. 1\ that any 

limit point r ro i = {^] t }k>i £ fe >i C([0,T],C\) is such that, for every k > 1, 7^ > 0, and 7^ 
is symmetric w.r.t. permutations. In Proposition 16.31 we also show that 

Tr(l-A 1 )...(l-A fe ) 7 S<C fc (4.3) 

for every t £ [0, T] and every k > 1. Note that, for finite N, the densities 7^ do not satisfy estimates 
such as (|4.3p (at least not uniformly in N), because they contain a short scale structure. Only after 
taking the weak limit, we can prove (|4.3|) . 

Step 2. Convergence to the infinite hierarchy. In Theorem 17.11 we prove that any limit point 
Foo.t = {l^]t}k>i £ 0fe>i c ([0,T],Cl) of F N;t = {7^1}fc>i with respect to the product topology 
T prod is a solution of the infinite hierarchy of integral equations (k = 1,2,...) 

7^= W (fc) (*)72; -»^oE f ^ k \t - s)Tr k+l U{ Xj -x k+l ) n ^] (4.4) 

3=1 Jo 
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with initial data 7^ = \<p){p\® k . Here Tr^+i denotes the partial trace over the (k + l)-th particle, 
and Li^ k \t) is the free evolution, whose action on /c-particle density matrices is given by 

W (*)( t ) 7 W := ettEj-iA^JfeJg-itEkiAi _ 

Note that (|4.4p is the (formal) limit of the iV-particle BBGKY hierarchy (|3.30p (written in integral 
form) if we replace the limit of NVn(x) with 8nao5(x) (see (|3.38p ). 

The one-particle wave function <p was introduced in (|2.14p . From (|2.13p and the positivity of the 
potential we note that 

CN > ty N , (H N + N) ^ N ) > NTr (1 - A) 7$ . (4.5) 

Since by (|2.14|) . 7^ — > \(p){(p\ as N — > 00, w.r.t. the weak * topology of £ 1 (L 2 (]R 3 )), it follows from 
(143]) that Tr (1 - A)|<p)(p| < C, and therefore that tp G i^QR 3 ). 
We remark here that the family of factorized densities, 

i k) =\n)^r, (4.6) 

is a solution of the infinite hierarchy (|4.4p if ^ is the solution of the nonlinear Gross-Pitaevskii 
equation (I2.16P with initial data cpt=o = tp. The nonlinear Schrodinger equation (|2.16p is well posed 
in i? 1 (R 3 ) and it conserves the energy, S(tp) := \ f |V</?| 2 + 47rao f \tp\ 4 . From tp G // 1 (R 3 ), we thus 
obtain that tp t G ^(R 3 ) for every iGK, with a uniformly bounded H -norm. Therefore 

Tr (1 - Ai) ... (1 - A k )\p t )(p t \® k < \\ n \\ k Hl < C k (4.7) 

for all t G R, and a constant C only depending on the i/ 1 -norm of tp. 

Step 3. Uniqueness of the solution to the infinite hierarchy. In Section 9 of [9] we proved the 
following theorem, which states the uniqueness of solution to the infinite hierarchy (|4.4p in the space 
of densities satisfying the a priori bound (|4.3p . The proof of this theorem is based on a diagrammatic 
expansion of the solution of (|4.4p . 

Theorem 4.1. [Theorem 9.1 of J3j] Suppose T = {j^}k>l £ 0fc>i C\ ^ su °h that 

Tr (1 - Ai)... (1 - A fe )7 (fe) < C k . (4.8) 

Then, for any fixed T > 0, there exists at most one solution T t = {7t }fc>l G C([0, T], 0/ 

\4-4\) such that 

IV-(l-Ai)...(l-A fc ) 7t (fc) <C* (4.9) 
/or a// i G [0, T] and /or all k > 1. 

5iep ^. Conclusion of the proof. From Step 2 and Step 3 it follows that the sequence T^,t = 
{T/vtH>i e © fe>1 C([0,T],4) is convergent with respect to the product topology r pro d; in fact 
a compact sequence with only one limit point is always convergent. Since the family of densities 
Tt = {7t }fc>i defined in (|4.6p satisfies (|4.7p and it is a solution of (|4.4p . it follows that r^t — > T t 
w.r.t. the topology r pro d. The estimates are uniform in t G [0,T], thus we can also conclude that 
%(7iVt> 7t ) ~^ 0- I n particular this implies that, for every fixed k > 1, and i G [0, T], 7^ — ► 7 t 
with respect to the weak* topology of Cl. This completes the proof of Theorem 12.11 Actually, the 
estimates are uniform in t G [0, T], and thus we can also conclude that rjkilNplt ; ) — ► 0. □ 
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Next we prove Theorem 12. 2\ to this end we regularize the initial wave function, and then we 
apply the same arguments as in the proof of Theorem 12.11 



Proof of TheoremlgM Fix k > and x 6 Cq° (R) , with < x < 1, x(s) 
x(s) = 0ifs>2. We define the regularized initial wave function 



1, for < s < 1, and 



X (kH n /N)i(j n 
\ X (kH n /N)iPn\ 



and we denote by ipN,t the solution of the Schrodinger equation (12. 2|) with initial data i/jn. Denote 
by ^N,t = {TjvtlfcLi the family of marginal densities associated with tpN,t- By convention, we set 

~(k) ' 

>N t := if A; > N. The tilde in the notation indicates the dependence on the cutoff parameter k. 
In Proposition l8.lt part i), we prove that 



k i\ T k 



{tl>N,t,H$r*l> N ,t)<C*N 



(4.10) 



if k > is sufficiently small (the constant C depends on k). Moreover, using the strong asymptotic 
factorization assumption (12.180 . we prove in part iii) of Proposition 18.11 that for every £ JCk, 



(4.11) 



as N — > oo. From (j4.10|) and (|4.1ip . we observe that the assumptions fl2. 13f> and (I2.14|) of Theorem 

~ ~(k) 

12.11 are satisfied by the regularized wave function ipN and by the regularized marginal densities *y N t - 

Therefore, applying Theorem 12.11 we obtain that, for every i £ 1 and k > 1, 



where <pt is the solution of (|2.16p . 



~(&) 
TNt 



(k) 

It remains to prove that the densities T N t associated with the original wave function ijjN,t (without 

~(k) 

cutoff k) converge and have the same limit as the regularized densities *y Nt - This follows from 
Proposition 18-H part ii), where we prove that 

UN,t-i>N, t \\ = Un-^n\\ <Ck 1 ' 2 , 

where the constant C is independent of iV and k. This implies that, for every £ we have 

^jW( 7 g-7g)|<^ 1/2 (4-13) 

where the constant C depends on j( k \ but is independent of N, k or k. Therefore, for fixed k > 1, 
t e R, J (k) e K k , we have 



(4.12) 



Tr J( fc ) (4\ 



< 



Tr 



(fe) 
TNt 



~(k) 
TNt 



< Ck 1 ' 2 + Tr J( fe ) 



~{k) 

TNt 



Tr J« (7$ 



\<Pt)(V>t 



\<Pt)(<fh 



(4.14) 



Since k > was arbitrary, it follows from (|4.12|) that the l.h.s. of (|4.14p converges to zero as N — » oo. 
This completes the proof of Theorem 12.21 □ 
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5 Energy Estimates 



In this section we prove two energy estimates that are the most important new tools used in the 
proof of the main theorem. Both estimates concern the smoothness of the solution -i/>jv,i( x ) of the 
Schrodinger equation (|2.2p . uniformly in N (for N large enough) and in t € R. However, due to 
the short scale structure of the interaction, Vjy, uniform smoothness, say in the x\ variable, cannot 
be expected near the collision points \x\ — xA ~ 1/N, j = 2,3, ...,N. The key observation is 
that x\ — * ^jvt( x ) wm nevertheless be smooth away from these regimes, whose total volume is 
negligible. For technical reasons, the excluded regime will be somewhat larger, \x\ — Xj\ > £, but 
still with Nl 3 <C 1. The same statement holds for the smoothness in an arbitrary but fixed number 
of variables, x\,... ,Xk- This is the content of our second energy estimate Proposition 15.31 

Our first energy estimate, Proposition 13. 1\ controls only two derivatives, but it is more refined: 
it establishes smoothness of V>7V,t( x ) in the X{ and Xj variables (for any fixed pair i,j) after removing 
the explicit short scale factor (1 — w(xi — Xj)). This factor represents the short scale effect of the 
two body interaction Vjv(zj — Xj) on the wave function and it is responsible for the emergence of the 
scattering length (|2.9|) . 

5.1 H"^ Energy Estimate 

In this section, we shall prove Proposition 13.11 We first collect some important properties of w(x) 
(12.121) in the following lemma. This lemma is an improved version of Lemma A. 2 from {8J. By 
defining p somewhat differently (see (|2.7|) ) . we also correct a minor error in (A. 6) and (A. 19) of [8j. 

Lemma 5.1. Suppose V > is smooth, spherical symmetric, compactly supported and with scattering 
length oq. Let 

poo 

p = sup r 2 V{r) + / drrV(r) (5-1) 
r>0 Jo 

and let a = a^/N be the scattering length of the rescaled potential Vn- Then the following hold with 
constants uniform in N. 

i) There exists a constant Cq > 0, which depends on the unsealed potential V , such that 

Cq<1-w{x)<1 for all x eM?. (5.2) 
Moreover, there exists a universal constant c such that 

1 - cp < 1 - w(x) < 1 for all x G R 3 . (5.3) 

ii) Let R be such that supp^ C {x £ R 3 : \x\ < R}. Then 

w(x) = j— r for all x with \x\ > R/N . 

Hi) There exist constants C\, C%, depending on V , such that 

\Vw(x)\<dN, \V 2 w{x)\ < C 2 N 2 , forallxGR 3 . (5.4) 
Moreover, there exists a universal constant c such that 

\Vw(x)\ < c-^tt, \Vw{x)\ < c-^-, \V 2 w{x)\ <c-^ for all x eM 3 . (5.5) 

T* If IT " 
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iv) We have 

87ra = J dxV N (x)(l-w(x)). (5.6) 

Proof. We prove part i) and iii) in Appendix [Dj Part ii) follows trivially by the definition of the 
scattering length o and by the fact that the potential has compact support. As for part iv), note 
that, due to the spherical symmetry of Vn and w(x), with the notation r = \x\, the function 
g(r) := rf(r) = r(l — w(r)) satisfies 

-9"(r) + \v N {r)g{r) = . 
By ii) of this lemma, g(r) = r — a for r > Ra. We thus obtain 

poo roo 

dx Vn(x)(1 — w(x)) = 4tt / dr r 2 Vjsr(r)(l — w(r)) = 8ir / drrg"(r) 

) (5.7) 



8tt lim (r g (r) — g(r)) |£ = 87ra . 



□ 



Proof of Proposition ^. 1\ For j = 1, . . . , N, we define 

ffc := -Aj + - V n(xj ~ xt). (5.8) 

Then we clearly have 

N 
3=1 

Since ip is symmetric with respect to permutations, we have 

N 

(if>, H%^) = Mrf) = N(N - 1)(^, f)lfoi>) + N(ji, f)fy) > N(N - 1)(^, di^) . (5.9) 

i,3 

Of course, instead of the indices 1, 2 we could have chosen any i ^ j. 
We have 

fhV = -AiV + ^V N (x 1 - x 2 )ip + - Y V n(xi ~ Xj)ijj (5.10) 

i>3 

Next we write ip = (1 — 1^12)^12 an d we observe that 

- Ai[(l - w 12 )<h.2] = (1 - wi2)(-Ai0i 2 ) + 2Vw 12 Vi0i2 + Aw 12 <Pi2- (5.11) 

Hence 

(1 - Wl2) _1 f)l [(1 - W 12 )(f>12] = ~ Ai^i2 + 2 VWl2 Vi0i2 

1 - t«12 

(-Ai + (l/2)Vjy(xi -x 2 )) (l-^ia) ^ , x 

+ ; 9i2 (5.12) 

1 — 1012 



2 Z_ — 1 ^7^12- 

i>3 
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Using the definition of w{x) (see (|2. 12|) ) . we obtain 

(1 - U>12) -1 I)1 [(1 - ^12)012] = ^1012 + 2 V ^( x l ~ x i)^i2 ( 5 - 13 ) 

i>3 

where we defined 

Li :=-A 1 + 2-^ Vl . 

1 — W\2 

Note that this operator is symmetric with respect to the measure (1 — wi 2 ) 2 dx, i.e. 

(1 - w 12 ) 2 (L lX ) = J (1 - w 12 ) 2 (Lx0) X = f (1 - ^i 2 ) 2 V^Vix • (5.14) 
Analogously to (|5.13p . we have 

(1 - W 12 )~ X §2 [(1 - Wl2)<^12] = ^2012 + 2" Vn ( X<1 ~ X i)<h2 ( 5 - 15 ) 

i>3 

with 

Vll?21 

L 2 = -A 2 + 2- — V 2 . 

1 - W12 

Therefore, from (15.91) we find 



'12 



N(N-1) J (I - w 12 ) 2 L 1( j) 12 L 2 cPi2 
N(N - 1) 



+ 



2 

i>3 



y^fjl- wn) 2 {V N (x 2 - xj) Li0 12 0i2 + Viv(a;i - Xj)4> 12 L 2 (j) 12 } 



N(N - 1) ^ /• . . 2 r . . . 

H 1 2^ I V l ~ w iV V N (xi - Xj)V N (x 2 - Xi)\<f>i 2 



2 



i,i>3 • 

N(N-1) I (1 - w 12 ) 2 L{0 12 L 2 cp 12 



+ X) ^ / (1 - ^i 2 ) 2 {Vjvfa - x^lV^I 2 + ^( X2 - Xj)|Vi0i2| 2 } 

i>3^ 

+ iV(iV-l)^ J (l-w 12 ) 2 V N (x 1 -x j )V N (x 2 



Xi)\4>12\ 2 

i,j>3 ' 

>N(N-l) J {I - w 12 ) 2 L{4> 12 L 2 cp 12 . 

(5.16) 

Here we used that the potential is positive and that the sum Ylj>3 Vn(xi — Xj) is independent of x 2 
(and analogously Ylj>3 Vn{x 2 — xj) is independent of x±). 

From (|5.16p we find 

(1>,H%il>)>N(N-l) J (l-^i 2 ) 2 Vi0 12 ViL 2 ^ 12 

= N(N-1) J (l-^i 2 ) 2 |ViV 2 12 \ 2 + N(N-l) J (l-^i2) 2 V 1 ^ 12 [Vi,L 2 ] ( /» 12 . 

(5.17) 
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To control the last term, we note that 



V 



Vw 2 l 
1' 1 

1 - w 2 i 



< 



|V 2 w; 2 i| / Vwu 



+ 



1-W\2 V 1_u, 12 

by ()5.3p and (15. 5p . for p small enough. Therefore we have 



< cp-. 



Fl ~ x 2| 



y (i-wi 2 ) 2 Vi0i 2 [Vi,l 2 ]0 12 <cpy 



< cp (l- w 12 y 



1 



< 



1 



Fl ~ x 2\ 
2 



|Vi^i2||V 2 9>12| 



Fl — X2\ 



Vl^l2 



(5.18) 



< C p y |ViV 2( />i 2 | z 

< cp ^ (l-w;i2) 2 |ViV 2 <; 
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where we used (|5.3p to remove and then reinsert the factor (1 — u>i 2 ) 2 (assuming p is small enough), 
and where we used the Hardy inequality to control the l/|x| 2 singularity. From (|5.17j) we have 



{ijj, H 2 N i>) > (1 - cp)N(N - 1) / (1 - w 12 f |ViV 2 </> 



This completes the proof of the Proposition 13.11 



>12\ 



(k) 

For fixed 2 < k < N and i,j < A;, with i ^ j, we define the densities l)q \j t by 

T&.t == (1 - ^i) _1 7^1(l - WijT 1 , 



(5.19) 
□ 

(5.20) 



where (1 — Wij) 1 = (1 — w(xi — Xj)) 1 is viewed as a multiplication operator. The kernel of 
is given by 



i,j,t 



t(*k; x'fc) = (1 - Ufa - Xj)) 1 (1 - w(x'i - Xj)) VjJ-tC**;**) • 



(5.21) 



Then, for every k, and every i,j < k, with i 7^ j, J^\j t is a positive operator, with Tr7^ • t < C, 
uniformly in JV, 

Proposition 5.2 (A-priori bounds for 7^ 3 - t )- ^or any sufficiently small p, there exists a constant 
C > 0, such that 

Tr(l-A i )(l-A j )^l >j>t <C (5.22) 
for all t 6l, 2 < k < N , i,j<k, i 7^ j, and /or a// ./V /arc/e enough. 

Proof. For fixed i / j we define the function 0ij,t by ^jv,t = (1 — w ij) ( Pi,j,t (the iV dependence of 
is omitted in the notation). Then we observe that 



Tr (1 - A,)(l - Aj) 7 g J)t = HSiSj fef = ||^,t 
with S n := (1 - An) 1 / 2 . Next we note that, by (jOj) . 

dx|^(x)| 2 < dx|^ t (x)| 2 < C 



2 + 2||V^ i .,i|| 2 + ||V i V^ 2 



i v jVi,j,t\ 



(5.23) 



(5.24) 
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uniformly in N and t. Moreover 

^iV,t(x) 



V 



i<Pi,j,t\ 



dx 



< / dx 



V.; 



1 — w(Xi — Xj) 

1 



ViVjv,t(x)| 2 + / dx 



Viw(xi - xj) 



(1 - w(Xi - Xj)) 2 



(5.25) 



(1 - w(Xi - Xj)) 2 

<C f dx|V^Ar )t (x)| 2 + C / dx- ( x 

J J \Xi-Xj\ Z 

<C J dx|V 4 V^(x)| 2 . 
where we used (15. 2p . (I5.5P and Hardy inequality. Next we note that, for every i = 1, . . . , N, 

{if) N>t , H N ij, Ntt ) > N(i/> Ntt , A^JV.t) = N J \Vii/jN,t\ 2 ■ (5.26) 

Therefore, from (|5.25ft . 

||V^|| 2 < CN-^iPn^HniPnj) = CN-\tI> n ,H n ,I) N ) < C (5.27) 



by (|2.13p and by conservation of energy. Finally, to bound the last term on the r.h.s. of (|5.23p . we 
note that, for a sufficiently small p, 



|V,;V 



i V j<Pi,j,t\ 



< C j dx (1 - w( Xi - Xj)) 2 iViVj&j.tCx) 



C 



N(N - 1) 
C 

N(N - 1) 



(5.28) 



{i> N ,H 2 N ^ N ) < C 



for all N large enough. Here we used (|5,2p in the first line, Proposition 13.11 in the second line, the 
conservation of in the third line, and the assumption (|2,13p in the last inequality. Proposition 
EJnow follows from (IQ3|) . l|533j> . (IQTj) . and (I5T28D . □ 

5.2 Higher Order Energy Estimates 

We will choose a cutoff length scale £. For technical reasons, we will have to work with exponentially 
decaying cutoff functions, so we set 

/i(x):=e z — . (5.29) 



Note that h ~ if Ixl S> A and /i ~ e if Ixl <C ^. For i = 1 



, A" we define the cutoff function 



■■= exp 



(4x>< 



Xj ) 



(5.30) 



for some e > 0. Note that 0j(x) is exponentially small if there is at least one other particle at distance 
of order £ from Xi, while #j(x) is exponentially close to 1 if there is no other particle near Xi (on the 
length scale £). 

As for the choice of £, to make sure that the presence of particles at distances smaller than £ 
from Xi is a rare event, we will need to assume A"^ 3 <C 1. This condition is not used in Proposition 
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15.31 below, but if N£ s 3> 1, then our estimates were empty in the limit N — » oo as the r.h.s. of the 
estimate (|5.33p below tended to zero. On the other hand, choosing £ too small makes the price to pay 
for localizing the kinetic energy on the length scale £ too high. In Proposition 15.31 we will actually 
have to assume N£ 2 3> 1. 
Next we define 

0< W (x) := = exp ^ J2 h & ~ x i)) ( 5 - 31 ) 

and their cumulative versions, for n,k £ N, 

0^ (x) := *W(x) . . . ^(x) = exp ( Yl E h fa ~ *i) | • ( 5 - 32 ) 

To cover all cases in one formula, we introduce the notation 0^ = 1 for any k < 0, n G Z. We will 

need to use the functions #| n ^ (instead of #i(x)) to take into account the deterioration of the kinetic 
energy localization estimates. For example the bound |V,-#j(x)| < C£~ 1 6i(x) is wrong, while 

|Vi^ n) (x)| < cr^-^Cx) 

(n) 

is correct and similar bounds hold for 0i . This, and other important properties of the function 
used throughout the proof of Proposition 15.31 are collected in Lemma lA.ll of the Appendix. 

Proposition 5.3 (H k energy estimates). Suppose £ 3> N~ 1 / 2 and that p (from (|2.7|) ) is small 
enough. Then for C$ > sufficiently small (depending on the constant (1 — cp) in Proposition \3.1]) 
and for every integer k > 1 there exists Nq = No(k, Co) such that 

(V, (H N + N) k $ > C k N k J <df\ | Vi . . . V fc Vf 

+ C%N k ~ l j ef^ I V?V 2 . . . V fe _ 1? /;| 2 (5.33) 

+ c%N k+1 j e^^x) v^a* - xk+i) I Vi . . . v fe _^(x)| 2 dx 

for every wave function if) S L 2 (\R. 3N ) and for every N > JVo- 

In order to keep the exposition of the main ideas as clear as possible, we defer the proof of this 
proposition, which is quite long and technical, to Section [9J at the end of the paper. 



6 Compactness of the Marginal Densities 

In this section we prove the compactness of the sequence rV,t = {7^1}fe>i w.r.t. the topology r pro d. 

(See Section H] for the definition of r pro d and recall the convention that 7^ = if k > N.) Moreover, 
in Proposition 16.31 we prove important a-priori bounds on ciny limit point i of the sequence Tnj- 

Theorem 6.1. Assume that p is small enough and fix an arbitrary T > 0. Suppose that Tff,t = 
{7jy^}fc>i is the family of marginal density associated with the solution ipN,t of the Schrddinger 
equation \2. 2\) . and that &2.13\) is satisfied. Then T^t G ©fc>i ^([^^li^ife) Moreover the sequence 
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Tnj G ® fe>1 C([0, T], is compact with respect to the product topology T pro d generated by the 
metrics rfk (defined in Section^. For any limit point Toot = { 7 oo\}fc>l> Too^t * s symmetric w.r.t. 



permutations, 7^\ > 0, and 



/or every k > 1. 



rr 7 S<l (6.1) 



Proof. By a standard choice of the diagonal subsequence -argument it is enough to prove the 
compactness of 7 )y-+, for fixed k > 1, with respect to the metric n^. In order to prove the compactness 

Ik) ^ (k) 

of 7]y t with respect to the metric we show the equicontinuity of 7 ]y t with respect to the metric r/^. 
The following lemma gives a useful criterium to prove the equicontinuity of a sequence in C([0, T],CV). 
Its proof is very similar to the proof of Lemma 9.2 in [5j; the only difference is that here we keep fc 
fixed and we consider sequences in C\, while in [8] we considered equicontinuity in the direct sum 
C([0, T],H) = ©fc>iC([0, T], Tik) over all k > 1, for some Sobolev space 

Lemma 6.2. Fix k G N and T > 0. A sequence ^ N t G N = fc, re + 1, . . 7)^ > and 
TV 7^ = 1 for all t G [0,T] and N > k, is equicontinuous in C({0,T], C^.) with respect to the metric 
r/fc, if and only if there exists a dense subset J k of Kk such that for any G J k and for every 
e > there exists a 5 > such that 



sup 

7V>1 



TrjW(l%! t -l$ s )\<e (6.2) 



/or a// 1, s G [0, T] witfi |i - s| < £. □ 



For the proof of the equicontinuity of "f N t with respect to the metric we will choose the set 
Jk in Lemma [6T2l to consist of all G /C^ such that SiSjJ^SiSj is bounded, for all i / j, and 
i,j < k. We recall the notation S n = (1 — A^) 1 / 2 . 

Rewriting the BBGKY hierarchy (|3.30p in integral form we obtain for any s < t 

7$ = 7$, -ijlf^r [-A, , 7 « ] -iJZtdr [V N (x t - Xj ), 7$ " 



:£/ dr[-A i)7 «]-i^ / dr^^-x,)/ 
i=l 75 i<j 75 

i(iV -k)J2 drTr k+1 [V N ( Xj - x k+1 ),^ 1) ] 

_• 1 ./s 



(6.3) 



Multiplying the last equation with G J7fc and taking the trace we get the bound (recall the 
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definition (|5.20p of the densities 7^ t ) 



Tr & (7$ " |<Er^ |Tr J^Sj - S^SJ^) SjJ&Sj 

k ,.t 

+ J2 dr Tr (SiSjjWSiSj) (s^S^VNixi - Xj )(l - w lj )Sr 1 Sj 1 
x (SiSjjflturSiSj) (.S; '.S' ; '(I ,r,/;.V ; '.S' ; •) 

ft „t 

dr|Tr (SiSjjWSiS^^SjHl-WiAS^Sj 1 



+ 1 



ft ~J 

£ y dr|Tr (57 1 ^ 1 W iV (x,-x fc+1 )(l-^, fc+1 )^ 1 ^ 



1 c-l 

j 



S k+1 Sji^ k+ljr SjS k+1 ) (Sj 1 S k ^ 1 (l - Wj jk +i)Sk+iS j 1 



+ (1- E jf H Tr (^ J(fc) ^) (sr 1 ^ 1 



1 c-l 



x (5 fc+ iS 7 -7^ti ) +ljr 5j5 fc+1 ) ^^JWjvCxj - x fe +i)(l - Wtf+JS^Sj x ) 



Here we used that Sfc+i commutes with J^ k \ Next we observe that (see Lemma 16.41 below) . 



\\Sr L Sr L NV N (xi - xj)(l - w^Sr'Sr'W < CN / V N (1 - w) < C, 
by part iv) of Lemma 15.11 Moreover 



\\S^Sj\l- Wij )S^Sj l \\<C 



(6.4) 



(6.5) 



(6.6) 



and 



\\ S j S k+l( 1 ~ w j,k+i) S k+i S j X H - S j S k+i( 1 ~ w j,k+i) S l+i S j 2 ( 1 - Wj,k+i)S k+1 Sj 1 



< C 



s j ^fe+iVfc+iCl - «; J -,fc + i)5 J . 2 (1 - u; J -, fc+ i)Vfc + i5^ 1 5 7 . 1 



+ 
<C* + 



^fc+^Vfc+i^fc+i)^- 2 (V k+ iw jt k+i)S k l 1 S j 1 



< c 



(6.7) 



In the last step we used the second bound in (|5.5p . Since G J7fe is such that H^Sj J^SiSj\\ < C 
for all i,j = 1, . . . , k, it follows from (|6.4p - (|6.7p that 



Tr ^ (fe) (7$ " 7^1) < C fc (t - a) max max sup Tr Stfj -y%\ , (6.8) 

V ' ' / n=fc,fc+l J<n r g[ s t] ' 
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for a constant depending on k and on j( k \ but independent of t, s, N. From Proposition 15.21 and 
from the fact that the subset J 1 ^ is dense in /C&, it follows that the sequence G C([0,T], 



is equicontinuous. Since, moreover, Tr J*rl = 1 uniformly in t G [0, T] and N, the compactness 

(k) ^ 

of the sequence j\t{ w.r.t. the metric rjk follows from the Arzela-Ascoli theorem. This proves the 
compactness of Tpf )t = {7^ t }fe>i G ©fc>i C([0, T], with respect to the product topology r pro d- 

Now suppose that Too,t = {l^\}k>i G ©fc>i C([0>^1>£fc) is a limit point of T^t with respect to 
r prod- Then, for any k > 1, 7^ t G C([0,T],£^,) is a limit point of The bound 



(fc) 

7oo,t 



< 1 



Tr 

follows because the norm can only drop in the weak limit 

To prove that 7^ is non-negative, we observe that, for an arbitrary ip G L 2 (]R 3fc ) with \\ip\\ = 1, 
the orthogonal projection P v = \(p)(<p\ is in ICk and therefore we have 

(ip, 7 £^) = Tr P^] t = lim Tr P, 7 gJ = lim (<p, <p) > , ( 6 .£ 



1 (AO (k) 

Similarly, the symmetry of 7^ t w.r.t. permutations is inherited from the symmetry of 7*^ for 



for an appropriate subsequence Nj with Nj — ► 00 as j — ► 00 

inherited frc 

finite N. For a permutation tt G 5fc, we denote by the operator on L 2 (M 3fc ) defined by 

E 7T (p(xi, ...x k ) = (p(x„i, ... , x nk ) ■ 

(k) 

Then the permutation symmetry of 7^, t is defined by 

^itW = lt]t (6-10) 

for every ir G To prove (|6.10p . we note that, for an arbitrary G /C& and a permutation 
7r G we have, for an appropriate subsequence Nj — > 00, as j — ► 00, 



Tr J« 7 g t = Hm J« 7 « = Hm Tr J<*>H, 7 «H^ = Km Tr H" 1 



(fc) 

. .. Nj.t 
TV- w-l r(fc)w ^( fc ) _ Tr ^( fe ) w-l 



where we used that, since G AC fc , also E~ l J^E n e IC k . □ 

In the next proposition we prove important a-priori bounds on the limit points r^t- These 
bounds are essential in the proof of the uniqueness of the solution to the infinite hierarchy (14. 41) . in 
Theorem 14.11 

Proposition 6.3. Suppose that p is small enough, and assume that \2. 13\) is satisfied. Let = 
{7oo\}fc>l £ ©fc>i C([0> T],£j,) is a limit point of the sequence T^t = {lNuk=l w - r -t- the product 
topology Tp r0( i. Then 7^ (has a version which) satisfies 

2V-(l-A 1 )...(l-A fc ) 7 £J t <C? (6.12) 
for a constant C\ independent of t G [0, T] and k>l. 
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Proof. We fix t as a function of N, such that N£ 2 3> 1, and N£ 3 <C 1- Moreover we fix e > so 
small that N£ 3 ~ £ <C 1. With this choice of £ and e, we construct, for integer n, k the cutoff functions 
Oj^(x) as in (|5.32p . For k E N, we will use the notation 

£> fc := Vi...V fc , := V'i . . . V' fe , with V$ = . 

We also set D k = I for k < to cover all cases in a single formula. From Proposition 15.3^ it follows 
that, for any fixed k > 1, 



I (6-13) 



C$iV fc 



for any N large enough (depending only on A;). In the last inequality we applied the assumption 

(k) 

For k = 1, . . . , N, we define the densities U N t by their kernels 

^l(x fc ;x' fc ) := y dxjv-fc 8jjf \x k ,x N _ k )Q k k \x' k ,x N _ k )D k i) N j(x k ,x N _ k ) D' k ip Nt (x k ,x N _ k ). 

(6.14) 

Note that the operator U^ t is the /c-particle marginal density associated with the iV-body wave 

function @ k k \]t)D k ipi\ij(x.). Therefore > 0. Moreover, it follows from (|6.13p that, for N large 
enough, 

Tr U nI = f H k) } 2 \D k ^N,t\ 2 < I ©fe-i \D k ^ N ,t\ 2 < C\ ■ (6.15) 

It follows from f)6. 15|) that for every fixed integer k > 1, and for every t £ [0, T], the sequence 

is compact w.r.t. the weak* topology of Ci. Moreover, if denotes an arbitrary limit point of 

Uft\, then 

Tr^S<C 2 fc . (6.16) 



Next we assume that 7^ £ C([0,T], C\) is a limit point of 7^ w.r.t. to the topology rj k . It 
follows that for any fixed t € [0, T], 7^ is a limit point of 7^ w.r.t. the weak* topology of C\. 
Because of the compactness of the sequence w.r.t. the weak * topology of C k , we can assume, 

(k) 1 (k) 

by passing to a common subsequence iVj, that there exists a limit point t £ C k of U N t such that, 

Tr J« 7 W -> Tr J<*) 7 <« (6.17) 

and 

Tr J« ^ - Tr tf^ (6.18) 

for every € K, k . For notational simplicity, we will drop the index i, but keep in mind that the 
limits hold only along a subsequence. 
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Next we fix G tCk such that Vi . . . VfcJ^V| . . . Vj is compact and such that 

sup / dx' fc £ l V < V 4l . . . Vi, . • • V; m J( fc )(x fe ; x' fc )| < oo 
Xfc 7 6=0 

r 4 

sup / dx fc ]T |V* n . . . ViXi • • • V; m J (fc) (x fc ;x' fc )| < oo 



for every j, m,n < k, and (ij, . . . ij), (ri, . . . , r m ) C {1,2,..., fc}. Then we have, applying (|6.17p to 
the derivatives of 

Tr Vi . . . V fe J^Vl ... VI jfil - Tr Vi . . . V fc j( fc V fc . . . V? 7 g t (6-20) 
as Ni — > oo. For such observable we rewrite the l.h.s. of (|6.18p . using (|6.14p . as 

TvJ^U^ t = /dx fc dx' fc dx w _ fc j( fc )(x fc ;x , fc )ef ) (x fc ,x JV -fc)©Af :) (4.XiV-fc) , , 

j (6.21) 

X -DfcV , Af,t(xfe,XAT_fc)£>fc'i/ , Ar,t(x / fc ,X Ar _ A .) . 

From (|6.2ip . we will show later that 

Tr J (fe) U$ t = J dx fc dx / fe dx A r_ fe (Ac D' k J (fc) ) (x fc ; x' fc ) ipN,t{ x k, x N-k)^N,t{^ Xjv_ fc ) + o(l) 

(6.22) 

as iV — > oo. 

Before proving (|6.22p . let us show how Proposition 16.31 follows from it. Equation (|6.22p implies 
that 

Tr J« U$ = Tr Vi . . . V* J<« VjJ . . . VJ 7 $ + o(l) 

(6.23) 

-^TrVi...V fc J«V^...VhS t 
as A?" — ► oo (using (|6.20p ). Comparing with (|6.18p . we obtain that 

Tr J<*> U« = Tr V X . . . V fe J« V* fe . . . Vj 7 £i • (6-24) 

Since the set of all G K, k with the property that Vi . . . V k J^V\ . . . V* k G K k and such that 
(|6.19p is satisfied is a dense subset of JC k , it follows that 

Vi • • • V kl £] t V* k . . . VJ = . (6.25) 

From (|6.16p . we find 

Tr (— Ai) . . . (— A fc ) 7( ^ t < Cf . (6.26) 

Now suppose that r Mi ( = { 7 <^t}fc>i £ C([0,T], C k ) is a limit point of the sequence Tjyj- Then, 

for every fixed k > 1 and i G [0,T], 7 £ t is a limit point of 7 jyl and thus satisfies (|6.26p . for a 
constant C2 independent of i and k. Moreover, for any m < k we also have 

Tr (- AO . . . (- A m ) 7 £ t < CJ». (6.27) 
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To prove the last equation, we repeat the same argument leading from (|6.14p to (|6.26|) . but with the 
densities replaced by 

^m,7V,t( X fc; X fc) = J dx^-fc & k k \x k ,X N _ k )®f\x k) X N _ k )D m 4> N j(x k ,X N ^ k ) D' m 1p N j(x k ,X N 



From (|6.26p . (|6.27j) . and from the permutation symmetry of 7^j\, we find 

Tr (1 - A X ) ... (1 - A,hS = E Q)tt (-A x ) . . . (-A m ) 7 g < (C 2 + l) k 



(6.28) 



(6.29) 



which completes the proof of Proposition 



It remains to prove (|6.22p . To this end, we rewrite the r.h.s. of (I6,2ip by using 0l fe ^ = 0^Q k __ l 
as follows: 

Tr J( fc ) U%1 = (I) - (II) (6.30) 

with 



(I) := J dxfedxfcdxjv-feJ^^Xfej-Xfe) Q^\(x k ,x N - k )G^ (x' k ,x N - k ) 

(II) ■= J dx fc dx' fc dx A r_ fc J (fe) (x fc ;x / fc )(l - 6»f ) (x fe ,x Ar _ fc ))G^ 1 (x fc ,x iV _ fc )e^ ) (x / fc ,x A r_ fc ) 

x D k ^ N)t (x k ,x N - k )D' k -ip Njt (x' k ,x N - k ) . 

By integration by parts 



(6.31) 



(/) = (la) + (lb) 



(6.32) 



with 



(/a) := - J dx fc dx^dx J v-feV fc J (fe) (x fe ;x^)e^ 1 (x fe ,x iV _fc)e^' ) (x / fc ,x i v_fc) 

X Dk-lTpN,t(*k,XN-k) D'k^N,t(^k^N-k) 

(lb) ■=- J dx k dx k dx N _ k J^ k \x k ]x' k )V k Q k k } 1 (x k ,x N ^ k )Q k k \x' k ,x N ^ k ) 

x Dk-iipN,t(xk,*N-k) D k ip N j(x' k ,x N _ k ) 
The main term is (la). To bound the term (lb), we use Schwarz inequality with some a > 0: 

|(I6)| < a J dx k dx' k dx N _ k \j( k ' ) (x k ;x' k )\ Vfee^ 1 (x fc) x iV -fc) \D k ^ Ntt (x k , x^^)] 2 

+ a" 1 J dx fc dx' fc dxiv-fc| J^(x k ;x' k )\ Q k k+1 \x k ,x N - k )\D k Tp Ntt (x k ,XN-k)\ 2 

< a (sup J dx' k \J^(x k ;x' k )\^ J dx \v kQ^l^ \D k ^ N , t (x)\ 2 

+ a~ l (sup J dx fc | J (fe) (x fc ;x' fc )| j J dx k dx N _ k ©^(x'fc, x N ^ k )\D k tl; Ntt (x k , xjv-fc 



(6.33) 



(6.34) 
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Using that 



we obtain that 



(6.35) 



dx 



(fc ^ \ ^ 

^ £ - x m )j ef+ 1) (x)p ifc _ 1 Vjv 1 t(x)| 

m=2 



(6.36) 



Xi - x m ) e^ + 1 1) ( x )lA-lV'Ar lt (x)p 



where we used the symmetry of the Dk-iipN,t w.r.t. permutations of the last N — k variables. Since 



~>k ™ 



E ( 7F E ^ " *») ) W < I^EE ^ - *m) I eE 1J (x) < Cefi(x) (6.37) 



(k+l), 



»( fc ) 



j>fc \ m=2 



i>fc m=2 



(see part ii) of Lemma lA.lj) , it follows from (|6.36p that 



Vk& [ k-i ' \ D k-ii>N,t\ 2 < cr 2 (iv - A;)- 1 / 6^ \D k ^ N)t 



»(*) 



< cr^jv - A:)- 1 J e k k _ 2 1] \D k ^ N>t \ 

< C k £- 2 (N - k)~ l 



(6.38) 



by (I6.13P (here the constant C k depends on k and on the observable J^). From (16.34p . from the 
assumptions (|6.19p . and again using (|6.13p . it follows that 



< C k (a{N - k)~ x r 2 + a" 1 ) = o(l) 

because N£ 2 > 1. 

Next we consider the term (II) in (16,311) , By Schwarz inequality, we have 
< a J dx fc dx' fc dx^ fe |j( fe) (x,.;Xfc)|0^^ 

+ a~ l J dx fc dx' fc dx^_ fc |j( fc )(x fe ;x' fe )|(l -^ fc) (x fc ,x JV _ fc ))Gf +1) (x' fc ,x w _^ 

X \D k lj}N,t{ x 'k,*N-k)\ 2 
<a(sup J dx' fc |j( fe )(x fc ;x' fc )|) J dx @^\x)\D k ^ N ^)\ 2 

+ a- 1 ( sup /dx fe |jW(x fc ;x , fe )|(l-^ fe) (x (fc ,x A r_ fe )) J 

x J dx'^dxjv-fe ® k k+1 \x k ,x N - k )\D' k 'il) Nt t(x' k ,XN-k)\ 2 

<C k [a + a- 1 sup /dx fc |jW(x fc ;x , fc )|(l-^ ) (x fc ,x 7V - fc ))) , 



(6.39) 



(6.40) 
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where we used (|6.13p . Next we note that 

! dx fc |J«(x fe ;x / fc )|(l-0f ) (x fe ,x JV _ fc )) < 2 —Y. \ ^k\J {k) {^ k )\h{x k -x m ) 

< C k Nl 3 -z f dx fe |V 4 fc j( fe )(x fc ;x' fc )| + |J (fc) (x fc ;x' fc )| 
because, with h(x) = exp(— (x 2 + 1 2 ) 1 / 2 /£), we have, by the Sobolev inequality, 

4 

I' dxh(x)\f(x)\ < ^Hill/Hoc < C£ 3 [J2 |V 6 /I • (6-42) 

J •* b=0 

From ()6.40p . (|6.41[) . and from the assumptions (16.190 we find 

\{H)\ <C k (a + a^iW 3 " 6 ) -> (6.43) 

as N —¥ oo, because N£ 3 ~ £ <C 1. 

From (|6.30p . (16.39P and last equation we find 

Tr J« U { *\ = ( dx fc d4dx^_ fc V fc jW(x fc ;x' fc ) ej^x*, x w _,) ef (x'^x^) 

7 (6.44) 

X D k -i1p Njt (x.k,XN-k) D'k^N^k^N-k) + o(l) 
Repeating the same arguments to move the derivative Vl from ^jvt to we obtain 



TrJ«t4 fe ;= /dxfedx'fedx^fcVfcVU^^x^x^efiCx^XAr-^eS^x^XAr-feJ 

j (6.45) 

X -Dfc-lV ; iV,t(xfc,X A r_ fe ) J D^_ 1 V'Ar it (x / A .,X A r_ A .) + o(l) 

Iterating this argument k—1 more times to move all derivatives to the observable, we prove (|6.22p . □ 

The following lemma was used in the proof of Theorem 16.11 and will also be used in the next 
sections, in order to bound potentials by the action of derivatives. 

Lemma 6.4. i) Suppose V G L 3 / 2 (R 3 ). Then 

dxV(x)\ip(x)\ 2 < C\\V\\ L3/ 2 I dx(\V<p(x)\ 2 + \ip(x)\ 2 ) (6.46) 



ii) Suppose V G L 1 (M 3 ). Then the operator V(x\ — X2), viewed as a multiplication operator on 
L 2 (M 3 x M 3 ,dxi dxa), satisfies the following operator inequalities 



V(x 1 -x 2 ) <C||F|| L i(l- Ai)(l- A 2 ), and V{x x - x 2 ) < C\\V\\ L i (1 - A x ) 2 . (6.47) 

The proof of (|6.46p is given in Lemma 5.2 of [8j, the proof of the first inequality of (|6.4Tj) is found 
in Lemma 5.3 of [10J. The last inequality follows from the usual Sobolev imbedding. □ 
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7 Convergence to the infinite hierarchy 

The aim of this section is to prove that any limit point r^t G fc>1 C([O,T],£j.) of the sequence 
T^t satisfies the infinite hierarchy 



Theorem 7.1. Suppose the assumptions of Theorem \2.1\ are satisfied and fix T > 0. Suppose 
Too,* = {l^]t}k>i G fe >! C{[Q,T],C\) is a limit point of Tnj = {lN,t}k=i with respect to the 
topology T pr0 (i. Then is a solution of the infinite BBGKY hierarchy 



7 2i = W«(t) 7 £J, - STriao^ / d.W«(f - s)Tr k+l \s( Xj - x k+1 ),^] (7.1) 

with initial data = I^X^I®* 1 . 

Remark. Note that in terms of kernels 

Tr fc+1 5(xj - x fc+ i)7^+ 1) ) (x fc ; x' fc ) = 7 ^+ 1) (x fe , x^-; x' fc , xj) . 



To define this kernel properly, we choose a function g G Cg°(K 3 ), g > 0, f g = 1, and we let 
g r (x) = r~^g(x/r). Then the definition is given by the limit 



= • 7co^s ( x fc ) 3y ) X fc ) ^ 



(7.2) 



The existence of this limit in a weak sense (tested against a sufficiently smooth observable) follows 
from the apriori estimate (|6.12p and from the following lemma (whose proof was given in Lemma 8.2 
in El). 

Lemma 7.2. Suppose that 8 a (x) is a function satisfying < 8 a (x) < Ca~ 3 l(|x| < a) and 
J 5 a (x)dx = 1 (for example 5 a (x) = a~ 3 g(x/a), for a bounded probability density g(x) supported 
in {x : \x\ < 1} ). Moreover, for G K k) and for j = 1, . . . , k, we define the norm 

I Ji% := sup (X!) 4 . . . <x fc )Vi> 4 • • • <4> 4 (| J (fc) (x fc ;x' fc )| + \V X] jW(x fc ; 4)| + \V X , J^(x k ; x' fc )|) 
x *> x fc 

, (7 - 3) 

for any j < k and for any function j( k \x k ; x' fc ) (Tie re (x) 2 := 1+x 2 ). Then if ^ k+1 >(K k+ i;x.' k+1 ) is 
the kernel of a density matrix on L 2 (R 3 ( k+1 ^), we have, for any j < k, 

J dx fc+1 dx' fc+1 J (fc) (x fc ;x' fc ) (8 ai (x k+1 - x k+1 )5 a2 {xj - x k+1 ) - 8(x k+1 - x' k+1 )5(xj - x k+1 )) 

x 7 ( fc+1 )(x fc+1 ;x' fc+1 ) 
< (const.)* p<-% (ai + v 7 ^) Tr {SjS^^ S 3 S k+1 \ . (7.4) 
Recall that S t = (1 - A^) 1 / 2 . TTie same bound holds if Xj is replaced with x'j in {1.1$ by symmetry. 
Proof of Theorem 7.1. For every integer k > 1, and every G IC k , we have 

sup TrjW ( 7 W t - 7 W)-,0 (7.5) 
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along a subsequence iVj — > oo. For an arbitrary integer A; > 1, we define 
In the following we assume that the observable E KL k is such that 



HS 



< oo, 



(7.6) 



where ||-A||hs denotes the Hilbert-Schmidt norm of the operator A, that is ||^4||hs = Trj4*A Note 
that the set of observables satisfying the condition (|7.6p is a dense subset of JC k . 
It is straightforward to check that 



\s 1 ...s k jWs 1 ...s k \\< nlj^n] 



HS 



Moreover, for any j < k 



< (const.) fc np^nl 



HS 



(7.7) 
(7.8) 



where the norm |||.|||j is defined in (|T.3[) . This follows from the standard Sobolev inequality 
(const.) ||/||^2,2 in three dimensions applied to each variable separately in the form 



< 



(^sup (x) i (x'r\V x J{x,x')\) < (const.) / dxdx' 



1 - A, 



(x) 4 (V x J(x,x'))(x') 4 



< (const.) Tr (1 - A)(x) 4 V J (x) 8 J* V* (x) 4 (1 - A) 

< (const.) Tr tfjn^TTl 7 



for 



with 0, = (x) + (1 — A) 1 / 2 . Similar estimates are valid for each term in the definition of 
j < k. Here we commuted derivatives and the weights (x); the commutators can be estimated using 
Schwarz inequalities. 

For G IC k satisfying (|7.6p . we prove that 



Trj( k )^ = TrJ^\ip)(ipf k 



and that, for t G [0, T], 



Tr J( fe ) 7 « = Tr J^ k \t)^\ - 8 W ^ / dsTr J^ k \t - s) 5( Xj - x k+1 ), ~ - 



i=i 



(fc+i) 



(71 



. (7.10) 



Note that the trace in the last term of (|7.10p is over k + 1 variables. The theorem then follows from 
(|7.9p and (|7.10p . because the set of 6 K. k satisfying (|7.6p is dense in K, k . 
The relation (|7.9|) follows from the assumption (|2.14|) and (|7.5|) . 

In order to prove (|7.10p . we fix t £ [0, T], we rewrite the BBGKY hierarchy (|3.30p in integral 
form and we test it against the observable JW. We obtain 



k . 

Tr J( fc ) 7 ^ =Tr (t)^ - if^ [ ds Tr J« W<*>(t - s)[V N ( 

i<j ^° 

-i(N-k)J2[ dsTrJ^U (k \t-s)[V N ( 
,=i ^ 



3)i 



(7.11) 



\ (fe+i)i 

Xj -X k+1 ),f N e J 
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From (|7.5p it follows immediately that 

Tr J( fc ) 7 « ^Tr J( fc ) 7( 

and also that 
Tr J^U^(t)4] Q = Tr (z/«(-t)J«) 7 « 



(fc) 

OQ.t 



(7.12) 



Tr (Ut k \-t)j^) 1 % = ^ J (k) U {k \t)Z\ (7-13) 



as N -> oo. Here we used that, if j( fc ) G £ fc , then also U^(-t)J^ G 

Next we consider the second term on the r.h.s. of (|7.1ip and we prove that it converges to zero, 
as N — » oo. To this end, we recall the definition (|5.2ip 



for every i ^ j, i,j < A;. Then we obtain 

Tr U^it- s)[V N { Xi - Xj),^. 



< 



Tr (SiSjiU^is - t)jW)SiS^ (^S7 1 ST 1 V N (x l - Xj )(l - w^S^Sj 1 ^ 



+ Tr (SiSjfldW^-QJ^SiS^ (sr^Sr^l-Wi^Sr 1 ^ 1 

x (SiSjy^SiSj) [s^Sj'v^-x^l-w^S-'ST^ 

Since, by part iv) of Lemma 15. 11 

\\Sr 1 Sr 1 V N (x l - Xj )(l - w lj )Sr 1 Sr 1 \\ <C J dxV N (x)(l - w(x)) < CN' 1 

and 

WS^sr^l-wixi-x^S^sr^Kl 

we find 

Tr J(*> U^it-s^VNixi-xj),^] < CN-^SiSj 



-i-l c-l 



(7.14) 



(7.15) 
(7.16) 



-t^S^WTrSfS]^. 

From \\SiSj (UW(S - t)jW) SiSj\\ = \\SiSjjM S l S j \\ < oo, and from Proposition [5J2] it follows 
immediately that, for any t £ [0, T], 



k nt 

V / dsTr J^U (k \t- s)[V N { 



Xi-Xj),^] -+0 



(7.17) 



as N — ► oo (the convergence is not uniform in k). 

Finally we consider the last term on the r.h.s. of (|7.1ip . First of all, we note that 



kf" [ dsTrJ^ k \t- s)[V N ( Xj - x k+1 ),^} 



(7.18) 
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as N — > oo. In fact 

Tr J«W«(i - s)\y N ( Xj - x k+1 ),^ 1] 
< 



Tr (SjUW(8 - t)J^Sj) (s^S^V^xj - x k+1 )(l - w jM1 )Sj£ 1 Sf 1 ') 
Tr (SjjWsA (sfSk+^l-wj^S-^S- 



.-i-l rr— 1 



-1 c-1 



(7.19) 



As in (I7TT5D we have [|S , J" 1 SfcJ: 1 Viv(xj - a?k+i)(l - wj^+^Sj 1 S^W < CN~ l . Moreover (see (fOl ). 

HSr^fc+^l - ^.fc+i)^^ 1 !! < C. (7.20) 
By an argument very similar to (|7.14p - (|7.17p and by Proposition 15.21 we obtain (|7.18|) . 



It remains to consider 

k ft 



ArV I dsTrjWuW(t-s)[V N ( Xj -x k+1 ) n % 
i=i J ° 



(k+lh 
s J 



k ,t 



3=1 



dsTr (U^(s - i)j( fc >) \NV N ( Xj - x k+1 )(l - w jtk+1 ),y^ 



j,k+l,s 



k ,.i 

-J2 ds Tr (^ (fc) (* " t)J (h) ) NV N (x 3 - x fc+ i)(l - w j!k+1 )^ k \ hs w j!k+1 
j=l J o 

k ,.f 

+ dsTr (^ (fc) (s-i 

.7=1 ^° 



(7.21) 



, i,*+i7Vj,fc+i,«( 1 ~ w j,k+i)NV N (xj - x k+ i) . 



The terms on the third and fourth lines converge to zero, as N — > oo. For example, the contributions 
on the third line can be bounded by 



Tr (U^(s - NV N (x 3 - x k+1 )(l - 

< \\Sj (U^(s - t) jW) ^WlSr 1 ^ (NV N ( Xj - x k+l )(l - w jtk+1 )) s^s^i 
x WS^Sk^wj^S^S^W Tr S]Sl +1 7 JJ \tJ +ljS . 
Then we use 



(7.22) 



(7.23) 



and 



||5j- Sk+iWj^+iSj S k+1 \\ < S k+1 Sj 1 Wj t k+iSl +1 w jt k+iS j S k+1 



1/2 



< ll^ 1 -^!^ 1 !! 172 + WS^SrWw^S^S^ (7-24) 
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To prove (|7.24p . we applied Lemma [6.41 and the fact that, by Lemma [5. 11 with R such that supp^ C 
{i£R 3 : |x| < R}, 



w(x) < C X (\x\ < R/N) + a X(|rE| > R I N ^ < C 



N\x\' 



and 



(the last bound is obtained interpolating the first bound in (|5.4p and the second bound in (|5.5p ). It 
follows that 



k ft 



s - t 



NV N (Xj - X k+ i)(l - Wj,k+l)lN^k+l,s W hk+l 



(7.25) 



<CtkN 1/4 max sup Tr SjS k+1 ^^V +1 SjS k+1 
se[0,t] 



which converges to zero, as N — ► oo, by using Proposition 15.21 The fourth line of (|7.21|) can be 
handled analogously. Hence, from (|7.2ip . 



NY] [ daTTj( k W e \t-a)\y N (x j 
i=i Jo 



x k +i),T N8 '\ 



/0 



j'=i 

+ C k T Ojv(l) 



(7.26) 



where ojv(1) as TV — > oo and C^y is a constant depending on k and on T. 

To handle the r.h.s. of (|7.26p . we choose a compactly supported positive function h £ Cq° 
with J dxh{x) = 1. For /? > 0, we define 8p(x) = P~ 3 h(x/(3), i.e. 8p is an approximate delta-function 
on the scale (5. Then we have 



o 



j'=i 



/ dsTr - t)J (fc) ) [Wtffo - a*+i)(l - u; iifc+ i) - 8vra ^(^ - Zfc+O.T/K- 



k+1,8 







+ J2 I dsTr (U^(s-t)j^) SnaoS^xj - x k+1 ),^ ] 



Y / dsTr {U^(s-t)J^) [87va S^ Xj - x k+1 ), 7 Jg 

■ i •/ 



ffc+1) 



+ C feiT (0(iV- 1 /2) + 0(^/2) 



(7.27) 



for some constant C k x which depends on k > 1, on T, and on J^-* (O^ 1 / 2 ) is independent of N). 
Here we used that, by (|5.6p . 

y dxWjv(x)(l - to(a;)) = 8vra , (7.28) 
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and we applied Lemma [7.21 To apply Lemma 17.21 we used Proposition 15.21 and that, by 

\tujf\s - t)j(% <Cpl 4 k \s - t)J<» Ql|| HS 

with a /c-dependent constant C. Since e~ l ^ s ~ t '^ (^"V^ - *-^ = (xj + 2(s — t)pj) m , for any j 
1, . . . , k, m E N, we obtain that 



imS fc) (* - t)f>% < c(i + |t - s\ 7 )\\nl j« nJH 



IIS 



To control the first term on the r.h.s. of (I7.27h we go back to Jpit ■ We write 



(h+i) 



~v"-r^ _ (*+!), ; i j -~y _i r_ 

'N,j,k+l,s ~ IN.s ' \ -i „,. M W,s ^ i 

\ 1 — U>j fc+l / ' 1 — Wj 



ffc+1) 



-7K 1) ( r - J ^ 



(7.29) 



When we insert (|7.29p in the r.h.s. of (|7.27p . the contributions arising from the last two terms in 
(|7.29p converge to zero, as N — » oo, for any fixed (3 > 0. For example, to bound the contribution of 
the second term on the r.h.s. of (|7.29|) . we use that 



Tr (uW(s - *)JW) toraoSpfa - x k+1 ), 



Wj,k+1 



1 ) 7E" 



< c 



+ CTr [U {k) (s-t 
-l 



5, 



fc+i 



1 - Wj,Jfc+i 



S k+i ) ( 5 JH-i7/£a {S^Mxj -x k +i)S k +i) 



fc+1 



ir o fc+1 7 Ar 



Now we have 



o_l ^j,fc+l c _i 

'-'fc+i-i „,, °fc+i 
1 — Wj,fc+i 



< CiV 



-i 



(7.30) 



(7.31) 



because w(x) < Ca\x\ 1 and thus, as an operator inequality, Wj^+i < CaS% +1 (and a ~ iV x ). 
Moreover 

^+17^ = {i>N,s, (1 - ^k+l)^N,s) 

< N-^n,,, (H N + iV)Vw, s > (7-32) 
= N- 1 (ip N ,(H N + N)?p N ) < C 

by the assumption (|2.13p . It is also easy to see that 

\\ S k+i S P( x 3 -Xk+i)Sk+i\\ < C(3~ A (7.33) 
for P < 1. The contribution arising from the last term on the r.h.s. of (17.291) can also be controlled 
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similarly. Therefore, it follows from (TT26j) . (TT27|) . (TT29|) . and (17301) that 



iVj^ / dsTrj( fc )^)(t- S ) V A r(x i -x fe+1 ), 7 j; 



(fc+i) 



o 



8vra ^ / dsTr (^( S -t)j( fc )j fyfo - x fe+1 ), 7, 



y (*+i) 

/oo,s 



(7.34) 



+ C k , T (o(p 1 / 2 ) + o N (l)) , 



where O7v(l)^0asA r ^oo (for any fixed (3 > 0). The first term is the main term. To control the 
second term, we rewrite it, for e > 0, as 







£/ d S Tr(wW( S -t)jW) ^(^-^O^^-T^ 



W,s '00, s 



- ]T f dsTrS^Xj - x k+1 ) 1 («< 1 >( S - t 
-^^d S Tr^-x fc+1 )(l- rT ^) 



_,(*+!) _ „(fc+l) 
W,s '00, s 



(fc+1) _ (fc+i) 
'iV,s '00, s 



'iV.s "00, s 



(7.35) 



The second term on the r.h.s. of (|7.35p can be bounded by using that 

Tr (U^\s - t)jW)6 ( Xj - Xk+1 ) (l - Y^s~) " * 

s - t)j( fc )) fyfo - s fc+1 )|| (Tr ^ +l7 ^+ 1) 5 fc+1 + Tr5 fe+l7 g+ 1 )S' fe+1 



< CfT 3 e (Tr S^gJ 1 * + TrSt +ll £ts 1] 

< Cp~ 3 e 



(7.36) 



where we used (17.32P and Proposition 16.31 Also the fourth term on the r.h.s. of (|7.35p can be 
controlled analogously. As for the first and third term on the r.h.s. of (|7.35p . we note that for every 
fixed e > 0, j3 > and s G [0,t], the integrand converges to zero, as N — > 00, by (17.51) . and because 

_t)j(*)) G /C fe+1 . 

(7.37) 



s - t)J ( > 8p{xj - x k+1 ) — — , 6p(xj - x k+ i)- 

/ 1 + eb k+ i 1 



+ eS, 



k+l 
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Since, moreover, the integrand is bounded uniformly in s S [0, t] (because for fixed e, (3 > the norm 
of the operators (|7.37p is bounded uniformly in s), it follows from Lebesgue dominated convergence 
theorem and from (|7.34|) that 



(fc+i) 

oo,s 



(7.38) 



k . 

= J2 [ dsTi (lA^(s-t)J^) [SnaoSpixj - x k+1 ), % 
3=1 J 

+ C k , T (<3(/3 1/2 ) + r 3 0(e) + ojv(1) 

where the convergence OAr(l)^0as./V^oo depends on e and (5. By applying Lemma 17.21 again 
and by using that, by Proposition 16.31 

max sup Tr(l-A J )(l-A fe+1 ) 7 2r ) ^ C - 
3=l,...,k te[0,T] 

we can replace 6p(xj — xu+i) with 5(xj — Xk+i) in (|7.38[) at the expense of an error 0(/3~ 1//2 ). 

From (Tmj) . (T?TT2D . (|715|> . (fTTTTI) . (|718|> . and (OH]) with tffo-orfc+i) it follows, letting iV -» oo 



with fixed /? > and e > 0, that 



j'=i 



^ / dsTv (u^ k \s-t)J^) \87Ta 5(x 3 -x k+ i),^]+O(^ 2 ) + ^O(e) . 



Eq. (|7.10p now follows from the last equation letting first e —* and then /? — > 0. 



□ 



8 Regularization of the Initial Wave Function 

In this section we show how to regularize the initial wave function t/'at given in Theorem [2? 



Proposition 8.1. Suppose that (2.11) is satisfied. For k > we define 

X (kH n /N)^ n 



\x(nH N /N)ip N \ 



is a cutoff function such that < x < 1> x( s ) = 1 f or < s < 1 and x( s ) 



Here x € Cg°( 

s > 2. We denote by j)^' , for k = 1, . . . , N, the marginal densities associated with ip^. 



(8.1) 
for 



i) For every integer k > 1 we have 



(i/) N ,Hfri/> N ) < 



2 k N k 



(8.2) 



ii) We have 



N 
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Hi) Suppose, moreover, that the assumption \2.18\) is satisfied, that is, suppose that there exists 
ip E L 2 (R 3 ) and, for every N 6 N and k = 1, . . . , N, there exists £jf~ fe) € L 2 (M 3(iV_fc) ) with 
W^N II = 1 sucn that 

. (8.3) 



lim Wn-^V^-QI 

N— >oo 



Then, for k > small enough, and for every fixed k > 1 and J^> G we /iGroe 

lim Tr - W){if\® k ) = . (8.4) 



Proof. The proof of part i) and ii) is analogous to the proof of part i) and ii) of Proposition 5.1 
in [SJ. Introduce the shorthand notation 3 := x( k Hn /N). In order to prove i), we note that 
l(-ffjV < 2N/k,)'B = 3, where l(s < A) is the characteristic function of [0, A]. Therefore 



(il) N ,H k N il) N ) 



fc 5^jy \ 

•^||' ^113^11/ 



< Hl^jv < 2JV/k)^|| < 



2 fc iY fc 



(8.5) 



To prove ii), we compute 

||3^7V - V^f = (l/>N, (1 - 3) 2 ^tv) < (V'iV, l(A6fTjv > iV)^Ar) • 

Next we use that l(s > 1) < s, for all s > 0. Therefore 

||H^jv - V'iv || 2 < —{^N,H N ip N ) < Ck 
by the assumption (|2.17p . Hence 

\\Zi/) N - iJ) N \\ < Ck 1 / 2 . 
Since HV'Jvll = 1, P ar t h) follows by (18. 8p . because 



(8.6) 



(8.7) 



JV 



3V 



v 



A' 



< \\ij) N - 3^tv|| + 

< 2||^jv-3^jv||. 



llp N 



3^ 



v 



\ZiPn\\ 



'} N - Eip N \\ + |1 



s^lll 



(8.9) 



Finally, we prove iii). For any sufficiently small n we will prove that for any fixed k > 1, € 
and e > (small enough) 



TrJ« 



< e 



1.10) 



holds if iV > No(k,e) is large enough. To this end, we choose ip* G .ff 2 (R 3 ) with ||y*|| = 1, such that 
||y-V*|| < e/(32A;||jW||). Then we have 



32||j( fc )| 



ill 



Therefore 

3^ E iff" ® « 



13(^0^ 



< 



< 4||^-^*®e^- fc) | 



5.12) 
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for k > small enough (by (|8.8p and because [|S|| < 1). Hence 
Eip N 



< 4||^ - ^ ® $- k) || + 4ii^* ^~ k) - v f k ® eif I 



< 



6||J( fe )| 



1.13) 



for N large enough. Here we used (|8.11|) and the assumption f)8. 3|) . Next we define the Hamiltonian 

JV N 



H N := - ^ A 3 - + ^2 V N (xi - Xj) 

j=k+l k<i<j 



U4) 



Note that Hn acts only on the last N — k variables. We set S := x{ k Hn/N). Then, from (|8.13|) . we 
will obtain 



7 4>N\ 



< 



3||J(*)|| 



J.15) 



for iV sufficiently large (if k > and £ > are small enough). 

Before proving (|8.15p . let us show how (|8,10p follows from it. Let 



v>f k ® $~ k) 



|H(^®4^ 



®k f -?jV 
V* ® ,,ft (AT-fe) , 

If? at 



since H acts only on the last N — k variables and since \\<p*\\ = 1. Moreover, we define 

7at (XfcJXj.) := / dx A r_ A .^ A r(x fc ,X A r_ A .)^ iv (x / A .,XAr_ fc ) . 



Note that is not symmetric in all variables, but it is symmetric in the first k and the last N — k 
variables. In particular, 7^ is a density matrix and clearly 



7« = i.e. 7 Ar fe) (x fe ;x' fc ) = n^(x J )^(4). 

Therefore, since \\ipN ~ •$n\\ < e/(3||J (fe) ||) by d8T5]) and since ||y>- < e/(32fc|| J (fc) ||), we have 

Tr (7? - i^r) I < |tt jw (jP - i^x^r) I + |tv (i^x^r* - i<^r 



1.16) 



< 2\\J^\\ \\i/j n -$ n \\ +2fc||J (fc) || || 



for N sufficiently large (for arbitrary k, £ > small enough). This proves (|8.10p . 

It remains to prove ()8. 15|) . To this end, we set ipN,* '■= ® > an< ^ we ex P aim the 

operator E — S = x( k Hn/N) — x( k Hn/N) using the Helffer-Sjostrand functional calculus (see, 
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for example, [3]). Let x De an almost analytic extension of the smooth function x °f order three 
(that is |c^x(- z )l ^ C| j/ 1 3 , for y = Imz near zero): for example we can take x(z = x + w) := 
\x(x)+iyx' \x)+X n '(x)(iy) 2 /2+x'" \x)(iyf /6]9(x,y), where 9 G C£°(M 2 ) and(9(x,y) = 1 for = x+iy 
in some complex neighborhood of the support of x- Then 



(3 -Sty 



AT,* 

7T 

Taking the norm we obtain 
Ck 



- I dxdyd 2 x(z) 



dx dyd z x{z) 



z - (kH n /N) z - ( K H N /N) 



AT,* 



1 



z-(kH n /N) 



(Hn — Hn) 



1 



z - (kH n /N) 



TV, 



(S - S)V>jv,*|| < 



iV 



dx dy 



\d*x( z )\ 



z-(kH n /N) 



[Hn — Hn) 



1 

z-(kH n /N) 



ipN,* 



Notice that the operator 



3.17) 



3.18) 



H N - H N = - ^2 A j + ^2 V N{x-i ~ xj) 

j=l i<k,i<j<N 

is positive hence (Hn — Hn) 1 / 2 exists. By using ||AB^>|| 2 < ll J 4|| 2 (V ; ) B*B\(j), we obtain 



z-(kH n /N) 



(Hn — H 



N 



z-(kH n /N) 



< 



(Hn — Hn) 1 / 2 - 



tpN,* 



1 



z -(kH n /N)\ 2 



(Hn - Hn) 1 / 2 



x ( ^ 



AT,*, 



(Hn — Hn)- 



-IpN,* 



3.19) 



3.20) 



z - (kHn/N) z - (kHn/N) 

Moreover (since H^yl 2 ^! = ||AB 2 A|| < ||v4C 2 v4|| for positive operators A, B,C with B 2 < C 2 ), 



(H N - H N ) 



1/2. 



\z-(kH n /N)\ 



: (H N - H N ) 



1/2 



1 



< 



< 



\z-(kH n /N)\ 
1 



(H N - H N )- 



z-(kH n /N)\ 



H 



1 



z - (kH n /N)\ " N )z - (kH n /N)\ 
CN 



\y\ 2 K 



3.21) 



for z in the support of Xi where we used the spectral theorem in the last step. On the other hand, 
the second factor on the r.h.s. of (|8.20p can be bounded by 



1pN,*: 



z-(kH n /N) 



(H n - H 



1 



N 



z-(kH n /N) 



< k(ipN,*, — — — — (-Ai + kV N (xi - x 2 ) + NV N (xi - x k+1 )) — L — — ip N ,*) ■ 

\ z-(kH n /N) z-(kH n /N) I 
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Here we used the fact that ipN* is symmetric w.r.t. permutations of the first k and the last N — k 
variables, and that the operator Hjy preserves this property. Since NVn(xi — Xfc+i) < C||V||£i(l — 
Ai) 2 , and kV N {xi - x 2 ) < C\\V\\ L i(l - Ai) 2 (see ([637])) we find 

1 / — 1 

ipN,*, — — — —{Hn - H N ) — — — - — -ip N . 

z — (kH n / N) z-(kH n /N) 

/ 1 1 \ ( ' 

* k (* N >* z-(«H N ,N) (" Al + (1 " Al) ^ z-(nH N /N) ^) ^I^M* 

because Ai commutes with H N (recall that ip N ^ = ipf k ® ^~ k) ). From (I8TT81) . (I8T20]) . (1831]) and 
(|HT22]> we find that 

||(E-E>jv,*|| ^C^N- 1 ! 2 

for a constant Ck )£ depending on k and e (through the norm ||<£>*||#2) but independent of k, for k 
small enough. This implies that 



< T" 7^1 ^7-7T^ IKS - S)^, 



s(^*®er fc) ) ii P^?*®^-*') ii " P(^*®er fc) ) ii (8.23) 



< 4 11(3-5)^*11 < 



6||J( fc )| 



for iV large enough (and assuming that e > and k > are small enough, independently of N). 
Here we used that (by (jPjl . ([8T8jh and (ISTTTI) ) 

us (Vf 8i e^- fc) ) ii > ii^ii - iis^ -iPn\\ - us - ^ ® e5T fc) ) II 

-ns(^®er fc) -^ fc ®er fc) )n (8-24) 



> 1 - Ck 1 ! 2 - o(l) > 1/2 

V ; 32||j( fc )|| ~ 1 

for K,e small enough and for N large enough. From (|8.23p and fj8. 13j) we obtain (|8.15p . This 
completes the proof of part hi). □ 

9 Proof of Proposition 15.31 

This section is devoted to the proof the Proposition 15.31 Let us recall the definition of the cutoff 
functions 

©i n) = ©i n) (x)=exp ( — ££*0*-^) 

from (|5.32p with the function h defined in (|5,29p . We introduce the notation hij = h(xi — Xj) and 
we also adopt the convention that ha = for any i E N. Moreover we recall that := Vi . . . . 

Proof of Proposition \5.3[ We prove (|5.33p by induction over k. For k = 1 we clearly have 

(iP, (H N + N)rp) > N I |Vi VI 2 + N{ - N ~ l) [ V N {x x - x 2 )|Vf • (9.1) 
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For k = 2 we have, from (|5.9p . (|5.16p (but keeping the term on the sixth line, which was neglected, 
because of its positivity, in the last inequality in (|5.16p ). (|5.17p . and (|5. 18|) we find, for p small enough 
(recall the definition of p in (|2.7|) ). 



>N{N-l){^M) + N^,^) 

> N(N - 1)(1 - cp) /(l - u; 12 ) 2 |ViV 2 0i 2 | 



(9.2) 



+ 



N(N — 1)(N — 2) 



(l-w 12 ) 2 V N (x 2 -x 3 )\V 1 <j )12 \ 2 + N / 



where fjj, for % = 1, . . . , N, was defined in (|5.8p . From the last term we get 

J \t)ii>\ 2 > J ( i ] \M\ 2 > J efAx^AiV + \Y.f ° { ? " X M + h -c-) (9-3) 



J>2' 



where h.c. denotes the hermitian conjugate. The last term is exponentially small in iV because on 
the support of the potential Vn{x\ — Xj) the point x\ is close to Xj (on the length scale TV -1 ) and 

this makes the factor 9\ exponentially small. Hence we find (with the notation V{ := <9^(j) where 

Xl = (Xi ,3?! ) G 



//i 3 /* 3 

^ 2) |Vi^l 2 + / E {(Vl^ 2) )(v0)viv{v + h.c.}+ / E viv{^ 2) viv;v{v 
ij=l ij=l 



-o(iV)<| / ^Wl 2 + / |Vl 



(9.4) 



by using |Vi#^| < CI 1 9[ l> from Lemma |A. 11 part iii). From part ii) and iv) of the same lemma 
we also have 



Vi0i 



(2) 



and therefore we obtain 

3 



2J = 1 



,00 



(Vl^ 2) )(v{v)viv{v 



< C£- 2 ^ 1} and V\df ) 



2/|(l) 



1 ) 



(9.5) 



( 2 )lV72„/,|2 . -1 / l^ 6 *! 



< (i) / ^ 2) iv 2 vi 2 + (iv) / ^ iJ ivi^r, 



(2) 1 2 



? f0 



-|ViVf 



(9.6) 



where we used that A^ 2 3> 1 (and an appropriate choice of the parameter a). Analogously 

3 



J IvivjtffVivvi^l <o(n) J ^ 1} |ViVf- 



(9.7) 
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From dOl-dO, we find 



(V, (H N + N) 2 tp) > N 2 (l -cp- o(l)) / (1 - ^i2) 2 |ViV 2( />i 2 



AT** 

+ — (1 - o(l)) / (1 - «;i 2 ) 2 V N (x 2 - x 3 )|V 1( /»i 2 | 2 (9.8) 



2 

" 3(2)|w2 ,,,2 at-2\ J / Z)(1)|V7 „/,|2 , //)(1)|„/,|2 



+ N(l-o(l)) J ef\V 2 ^\ 2 -o(N 2 )^J 0[ 1> \V 1 ^\ 2 + J o^M 

Next we apply Lemma 19.41 (with n = 0) to replace, in the first and second term on the r.h.s. of 
the last equation, cpi 2 by ip. We find 

(V, (H N + N) 2 ip) > N 2 (l -cp- o(l)) y ^ 2) |ViV 2 ^| 2 

+ ^(1 - o(l)) / 6f ] V N (x 2 - x 3 )|V^| 2 + N(l - o(l)) [ #f } |V 2 Vf (9.9) 



2 

-o(tv 2 ) y {^iv^P + ^IVI' + iv^^i-x,)^! 2 } . 

By (HH]) we have 

o(N 2 ) J ^ ) \V 1 ^\ 2 + e^ ) \^\ 2 + NV N (x 1 -x 2 M 2 } <o(N)(^,(H N + N)^) 

< o(l)(i>,(H N + N) 2 ^). 

Hence, from (19.91). we obtain 



(9.10) 



(9.12) 



(1 + o(l))(^, (H N + N) 2 ^) > N 2 (l -cp- o(l)) / ^IViVz^l' 

+ ^(1 - 1 - x 3 )|V^| 2 (9.11) 

+ N(l-o(l)) J 0$ 2) |Vfy| 2 . 
It follows that, for p small enough, there exists Co > such that we have 

(V, (H N + Nfi>) > C 2 N 2 J flf ^ViV^I 2 + C 2 N 3 J ef ] V N (x 2 - x 3 )|Vn/f 

if N is large enough. 

We assume now that fl5.33|) is correct for all k < n + 1 and we prove if for k = n + 2, assuming 
n > 1. To this end we note that, for N > iVo(n), using the induction hypothesis we have 

(H N + iV) n+2 ^} > (fl^, (H N + N) n H N ^) 

> c%N n J l^F^p (9 13) 

> c n iv ra y e(r+ 2 ) | A^Vf 
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where we used that 1 > 9^ > #| n+2 ) for every i = 1, . . . , re. We write Hn = Ylf=i tff^i w ^ n 

(n) = f -Aj + \ Ei >n ^j y N{xi - Xj) if j > n 

qj ' I -A,- + i E*<n VKrfa - Xj ) + E l>n Vn( Xi - Xj ) if j < n . K ■ > 

Then we have 

(4>,(H N + Nr +2 iP) >C r l N n Y, J ® { n +2) Dn^^D^f^ 

i,j>n 

+ CSN n lJ2 [Q { n +2) D^^Dn^^ + h.A (9.15) 

[ i<n<j J ) 
i,j<n 

The last term on the r.h.s. (where i,j < n) is positive and therefore it can be neglected. In the 
first term on the r.h.s. we can neglect all terms where i = j (because they are all positive). Therefore 
we obtain 



i,j>n,i^j ' 



(9.16) 



i<n<j 



In Proposition 19.11 below we give a lower bound for the first term in (|9.16p . while Proposition 19.5 
estimates the second term. Combining these two estimates, we find that, for p small enough (inde- 
pendently of N and n) and for N large enough, 



(V, (H N + iVT+V> > QiV n+2 (l - cp - o(l)) j ei n + 2) \D n+2 i;\ 2 

+ czN n+1 (i-o(i)) | ei n + 2) iVxAmVf 

+ (1 " O(l)) y - X„ +3 ) | A,+lVf 

where the error £l n (ip) is given by 

= o(N n+3 ) J G(T +1 ) V^v(x n+1 - x n+2 ) \D n Vf 

+ o(N n+1 )H e(r +1 )|v lJ D n </f + j e^jVi^-xVp 



(9.17) 



+ {n^)U e(r +1 )|zwi 2 + / ©iTilA^I 2 

+ / et^lDn-^l 2 + j eir 3 2) |^- 2 ^| 2 } 



(9.18) 
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Now we use the induction hypothesis, Eq. f|5.33j) . with k = n — l,n,n + 1 to bound the negative 
contributions. For example, (|5.33p with k = n + 1 implies that 



o( AP+3) J e (n+D VN{Xn+l _ Xn+2 ) \ Dn ^ < o(A f)(^ + AT+V) < o(l)^ s (#iv + N) n+2 ^) 
because Hn > 0. The other terms in (|9.18p are treated similarly. It follows that 

(1 + O(l)) (V, (H N + iV)™+V) > CJiV n+2 (l -Cp- o(l)) | | A^Vf 

+ C ^ +1 (l- O (l)) | 8i"+ 2) |ViA*iV| 2 

+ (1 " o(l)) J V N (x n+2 - x n+3 ) \D n+1 ^\ 2 . 

Thus, if p and Co are small enough (independently of n), we can find Nq(u + 2, Co) > No(n, Co) such 
that 

<V>, + AT +2 V>> > c " +2 iv^ 2 / ei"+ 2) |^ +2 ^| 2 + c "+ 2 Ar" +1 / ei n + 2) iVxiwvf 

(9.19) 

+ c n +27V n+3 /" Q J+ 2 )vi v (x w+2 -x n + 3 )|Z> M .i^ 2 . 



□ 

In the rest of this section we will state and prove Propositions I9.ll and 19.51 used in (19. 16H . Both 
proofs will be divided into several Lemmas. 

Similarly to the -ff^-energy estimate from Proposition [343 the key idea in Proposition 19. 1 1 is that 

(n) 

ip can be conveniently estimated by the derivatives of (pij, where <j>ij is given by the relation 

TP = (1 

— Wij)(pij. The estimates of all errors are done in terms of 4>ij and its derivatives. Finally, 
Lemma [9.41 will show how to go back from the estimates on </>y to estimates involving ip with a cutoff 
supported on a bigger set. 

Proposition 9.1. Suppose p is small enough and £ S> iV -1 / 2 . For i = 1, ... ,N, let r^™"* be defined 
as in \9. 11$ . Then 



i,j>n,i^j 



> c $N n+2 (i -cp- o(i)) j ei n + 2) \D n+2 ^< 



(9.20) 



+ — (i - o(i)) / ei n + + 2) v N (x n+2 - x n+3 ) \D n+ ^\ 2 - n n (ip) 

where the error term Q n (ip) has been defined in $9.18\) . 

Proof. For any i ^ j, i,j > n, we write ip = (1 — W{j)(f>ij. Then we have, similarly to (|5,13p . 



(1 - Wij) \H [(1 - Wij)(f>ij] = -Aifcj + 2 ^ Vifaj + ~ ^ Vzv(aJi - x mJ < , 



1-wu lTLJ 2 

(9-21) 



Li<t>ij + \ Y Vn ( 



m>n, m^i,j 
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where the differential operator Li := — Aj + is such that 

J(l- Wij f(L^) X = J{l- Wij )mL iX ) = j (1 - w^fV^ViX- (9-22) 
Note that the operator Lj also depends on the choice of the index j. Analogously, we have 
(1 - Wij)- 1 ^ [(1 - Wij^ij) = Lj^ij + - Yj ^M 3 



[Xi - Xm)(pij 



m>n, m^i,j 



with Lj = — Aj + 2 j^j". V j . Note that Z?„ commutes with Lj, Lj and 1 — Wij if i, j > n. The l.h.s 

m>n,my^i ,j 

\ L i + \ V N {xj-x r )\ 



of (|9.20p is thus given by 

C "iV" £ J (1 - G^ +2 ) £> n 



i,j>n,i^j 



(X m Xi) 



x £L 



> C^iV" £ y (1 - ^) 2 6^+ 2 ) Li-D„ <t> tj Ai v ,„ 



i,j>n,i^j 



+ 



C%N r 



J2 j{l-w ij fe^V N {x j -x r )L i D n 4> ij D n ct> ij ^h.c., 



i,j>n,i=£j r>n,r^i,j 



because of the positivity of the potential. Proposition 19.11 now follows from Lemma 19.21 and Lemma 

□ 



731 where we consider separately the two terms on the r.h.s. of the last equation. 
Lemma 9.2. Suppose the assumptions of Lemma \9.1\ are satisfied. Then we have 

C£N n j (1-Wij) 2 Li D n faj Lj D n cf>ij 

i,j>n,i^j 

> c£N n+2 (i -cp- o(i)) J ei n + + i 2) | A*rf 
-o(n^) (| eg* 1 * |zwi 2 + j e^ \ Dn ^' 

Proof. By the symmetry (|9.22p we have 



(9.23) 



i,j>n,i^j 



C%N n Y /( X " { n" +2) V^n ^ V l L J D n faj + V 4 ei n+2) V*D„ faj LjD n 



i,j>n,i^j 



C£N n /( X " w iS? { @ n +2) l V i V ^n ^l 2 + V i G^+ 2 ) ViD n ^j V 3 ViD n 4> i3 



+ Vi&l? +2) VjV t D n faj VjD n faj + V l V,e(T +2 ) ViD n faj VjD n fcj 

+ e^+v ViDnfaj T,./ v ;/;, f o (; } 



(9.24) 
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To bound the second and third term on the r.h.s. of (|9.24|) . we note that, by part hi) of Lemma lA.ll 



ve (n+2) 



< Ct 



1 E O e!," +21 

\ m=l / 



Therefore the second term on the r.h.s. of (|9.24p can be bounded by 



E 



i,j>n,i^j 



+ cr 2 a- 1 Y, I (I 



2n+2 



13 J 



m=l 



for some a > 0. Next we use that 0y = — w^) 1 . Since i,j > n, we have 



Vj-D n (-0(1 - Wij) x ) = ViWij(l - Wij) 2 D n ip + (1 - Wij) 1 VjD n V' 



and thus 



(9.25) 



(9.26) 



{l-WijflViD^A 2 < 2 



^pi-) |A^f + 2\VM\ 2 < - \D^\ 2 + 2|ViAWf • 

J- IV 4 <j / Xj X A 



(9.27) 



Therefore the second term on the r.h.s. of (|9.26|) is bounded by 



// on+2 n 
i, j >n,ij=j \ m=l 



Km ei,"+ 2 > |V,D„*j| 



^ E 

*° E 



in+2 



£fcj @i n+2) {|V^| 2 + r-^-T2lA 

m=l / 



e(™ +2 ) |v^| 2 + 



m=l 



(9.28) 



where we used Hardy inequality and the fact that i 7^ j and i > n. Using a bound similar to (|9.25p . 
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and part ii) of Lemma IA.ll we can continue this estimate 

r /nn+2 n \ 2 

£ J 0- - W ^ — £ k ^ Q n +2) I ViA, <t>i 
i,j>n,i^=j \ m=l / 



m=l 



m=l 



< 



G E/(^EE^) ©r 2) iv^ 

i>n \ j>nm=l J 



+ ct 



j>n m=l 
2 n+2 



EE^ 



2 n+l 



/ Z / Z / 

j>n m=l 



££M i n+2) i^p 



i>n m=l 
2 



< cE/ « n+1) l v ^Vf + cr 2 y ew |a^| 

i>ra 

because of the permutation symmetry of ip. From (|9.26p we find 



(9.29) 



t,j>n,i^j 



J{l- Wl] ) 2 VjQ^ +2) ViD n faj V J V l D n , 



< aN 2 J (I- w n+hn+2 ) 2 Q^ +2) \D n+2 <t> n+1>n+2 \ 

+ a - l cr 2 N I e(r +1 ) \D n+1 ^\ 2 + a~ 1 cr i I e[") \D n ^\ 2 



< o(N 2 ) (J (1 - w n+1 , n+2 ) 2 G(r +2 )|D n+2 </» n+1 , n+2 | 2 + J \D n ^\ 2 + j |A^I 

(9.30) 

for an appropriate choice of a (using that Nl 2 S> 1). In the last term we also used that 9^ < 1. 

The third term on the r.h.s. of (|9,24[) . being the hermitian conjugate of the second term can be 
bounded exactly in the same way. 

Now we consider the fourth term on the r.h.s. of (19.241) . To this end we use that, since i ^ j, 
and i,j > n, we have, by Lemma lA.ll part v), 



ViVj e^ +2) 



< CI 



-2 



' 2 n+2 n 

2^ h 



<«.< 

m=l 



' 2^+2 n 



71=1 / 



(9.31) 
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Therefore 



n Yi] 



r ( 9 n+2 n \ z 

i,j>n,iy^j \ m=l / 

< (iv 2 ) (| e(r +1 ) |zwi 2 + 1 e£>i lA^f 



(9.32) 



where in the second line we used (2.51) and a Schwarz inequality, in the third line we used the bound 
(f9T29|) . while in the last line we used Ni 2 > 1. 

Next we consider the last term on the r.h.s. of (|9.24p . To this end we note that, by (|5.3p and 



V- J — 



\V 2 w 



J'\ 



1 - Wi 



+ 



1 — w. 



1 3/^ 3-* j | 



assuming that p is small enough. Therefore, the terms in the sum on the last line of (|9.24[) can be 
bounded by using Hardy inequality as 



(1 - Wij ) 2 ei n+2 ) ViD n O . [V,, Lj]D n <fe 

<c P [(i-w^e^—L-^D 



< cp I T -—^\v l D n ^\ 2 



<c P / e^\VjViDM 2 + c 



nYi] I 



v, er 2) 



\ViD, 



(9.33) 



nY%3\ 



J97251 



< cp {I- Wij ) 2 |ViViD ft y | 



/• /on+l n \ 2 

+ ci- 2 hi- Wij f — Y, h jm & { r? +2) \ViD n < 

\ m=l / 



Next we sum over i, j > n (i ^ j); to control the contribution originating from the second term on 
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the r.h.s. of the last equation we use (|9.29p . We obtain 

2 | ji 1 ~ W ^ @ n +2) ^i D n [Vi, Lj]D n <P 

<c P Y, Jii-w^e^WjViD. 



i,j>n,i^j 



n ' 'ij I 



+ cr 2 N I e^ +1 ) | A+iVf + cr 4 / e£»> | a^I 2 (9-34) 



<C P ^ Jil-WijfG^lVjViDnfrjl 2 

i,j>n,i^j 

+ o(n 2 ) (| |A+iVf + f |AVI 

Inserting (IPO]) . and (I9341) into the right side of (19341) it follows that 

C£N n J (1 - ^) 2 e^+ 2 ) A y A ^ 

> C£N n+2 (l -cp- o(l)) y (1 - w nW2 ) 2 \D n+2 ^ n+1 , n+2 \ 2 (9-35) 

- o(N^ 2 ) ( [ | Am M 2 + [ ©£i I A,Vf N 



Lemma [9.2l now follows from (|9.41j) in Lemma f9.4l below that shows how to replace estimates involving 
the function fa = (1 — Wij)^ 1 -^ with estimates on tp. □ 



Lemma 9.3. Suppose the assumptions of Lemma \9.1\ are satisfied. Then we have 

Y j ( 1 - W V ) 2 « n+2) V " ( X 3 ~ X r ) L i D " $ij D " <t>ij + h.C. 



C%N n 



2 

i,j>n,i^j r>n, r^=i,j 

> ^^(1 - o(l)) J eiti 2) VV(x n+2 - x„ +3 ) \D n+1 ip\ 2 (9-36) 
- (iV" +3 ) / V N {x n+1 - x„,+ 2 )| AVf . 



C%N n 



Proof. Using ([<H2l) . we find 

^ y (1 - ^) 2 ei n+2 ) ^ (xj - s r ) UD n faD n fa 

- Vj r>n,r^i,j 

Y ^2 j - Wij) 2 V N (xj - x r ) 



2 

i,j>n,ij^j r>n,r^=i,j 

~ 2~ 

i,j>n,i^=j r>n, r^i,j 



(9.37) 
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Using (|9.25p (with j replaced by i) the second term in the curly bracket can be bounded by 



E E 

■,j>n,i^j r>n, r^=i,j 



J(l- WlJ ) 2 V N ( Xj - Xr) Vi®^ +2) V l D n hjDwy, 



N 

^ Ca E E I £ - Wo? & ( n +2) V N (Xj - X r ) IViD^ij] 



(9.38) 



W; 



J2h im ) e^V N ( Xj -x r )\D^ 



m=l 



Since i, j > n, and ip = (1 — u>ij)(f)ij, the second term can be estimated as 



// on+2 n 
Ve 

/ . ; 11 . I r i i • \i.r-pi. i V m=l 



him) & { n +2) V N ( Xj - x r ) |LWf 



<«-V^ £ /(^EE^) @^v N ( Xj -xr)\D n ^ 

j>nr>n,rj^j \ i>n m=l / 

< Ct~ 2 a x EE/ V N ( Xj - x r ) \D n ^\ 2 

j>n r>n, rjtj 

= Ct- 2 a~ l (N - n)(N - n - 1) / e^ +1 ) V N {x n+1 - x n+2 )\D n ^\ 2 , 



(9.39) 



because of the permutation symmetry of ip and Q^ +l \ From ()9.39p and (I9,38p . it follows that 

E E " W ^ V i @ n +2) V N (Xj ~ X r ) V l D n ^D n 4> 

" '*:3 

N 

<o(l) ]T E l0--mi) 2 G^ +2) V N (x j -x r )\V i D n <l> 



i,j>n,i^j r>n,r^i,j 



(9.40) 



+ o(N 3 ) [ @^V N (x n+1 -x n+2 )\D n iP\ 2 



where we used that N£ 2 3> 1 and we made a suitable choice of the parameter a. Inserting this 
bound into (|9.37p , using the permutation symmetry, and (|9.42p from Lemma 19.41 the lemma follows 
easily. □ 

The next lemma, showing how to replace estimates on fyj with estimates on ip, has already been 
used in the previous proofs. 

Lemma 9.4. Suppose the assumptions of Proposition [5731 are satisfied. Recall that (fiij is defined by 

ip = (l - Wij)4>ij. 



i) For n > 0, we have 

f (i-w n+1 , n+2 ) 2 ei n+ v \D n+2 <p n+hn+2 \ 



> (i - / < + i 2) I^WI 2 - o(i) { j e(" +1 ) |zwi 2 + / |A^<I 2 } • 



(9.41) 
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ii) For n > 0, we have 

J (1 - w n+ i jn+2 ) 2 6^ n+2) V N (x n+2 - \D n+1 4> n+ltn+2 \ 2 

> (1 - o(l)) y ei n + + 2) Fat(^ +2 - x n+3 ) |D n+1 ^| 2 - o(l) y ^(x n+1 - x n+2 ) \D n ^\ 2 . 

(9.42) 

Proof. In order to prove part i) we start by noticing that 

J (I- w n+1 , n+2 ) 2 @^ +2) \D n+2 <p n+hn+2 \ 2 > y (1 - w n+hn+2 ) 2 Gi n + + 2) \D n+2 <j) n+l , n+2 \ 2 . (9.43) 

Using that </> n+ i,n+2 = (1 - w n +i,n+2)~ 1 ip we find 

n a 1 n / 1 Vu, n+l,n+2 n „ , 

F> n +2<Pn+l,n+2 = Z Ai+2 W + Tx V? D n V n+2 V 

1 - «Wl,n+2 (1 - Wn+l,n+2) 

and thus, from f)5.3j) bounds] 

y(l - ^n+l,n +2 ) 2 O^ 2 * | J D n+2 n+1>n+2 | 2 

> y 8<f+ 2) | ZWVf -C j ®^ ] \Vw n+w \ \D n+2 ^\ \D n+1 xl>\ (9.44) 
-C J Q%+p (|V Wn+ i, n+2 | 2 + |V 2 ^ n+1 , n+2 |) \D n+2 xJ;\ \D n rP\ . 
The second term can be bounded by 
®n+l |Vw n+ i jn+2 | \D n+2 tp\ \D n+1 ip\ 

<aj Q { ^ 2) \D n+2 M 2 + a- 1 f Q^\Vw n+ i, n+2 \ 2 \D n+1 t/;\ 2 
< a f e££ 2) \D n+ ^\ 2 + a' 1 J ei n + + 2) x(Nn+l - x n+2 \ > I) |V™ n+1 , n+2 | 2 \D n+l ^\ 

+ / ©i+ + i 2) x(kn+i - x n+2 \ < e)\Vw n+1 , n+2 \ 2 \D n+1 ^\ 2 



2 

(9.45) 



fl("+ 2 )|n .„,/,|2 a. /-,v-l„2 /" o("+2) X(kn+1-Xn+2| 



> 



|3-n+l 2-n+2| 



+ Ca^iV 2 y 8^ X (K+i - x n+2 | < *) | An-iVf 

where in the last inequality we used that, by Lemma [5.1l |Vu; n +i )n + 2 | < CN . Moreover we used that 
Vw(x) = a/\x\ for \x\ > R/N (with R such that suppU c{i£l 3 : \x\ < R}), and that R/N < t 
for N large enough. Using that 

0%+? X {\x n+1 - x n+2 \ <£)< Ce~ cl ~ e (9.46) 



49 



we have (recall that a = ao/N) 
J ei n + + 2) |Vw n+ i, n+2 | \D n+2 iP\ \D n+1 il>\ 

Since N£ 2 > 1, we find 

ei" f + 1 2) |Vu; n+1 , n+2 | \D n+2 iP\\D n+1 ^\ < o(l) | j ei n + 2) \D n+2 ^\ 2 + J \D n+1 i>\ 2 \ . (9.47) 

As for the third term on the r.h.s. of (|9.44f) . we proceed as follows. 
J e { n+? (|Vw;n+i,n+ 2 | 2 + |V 2 u; n+ i, n+2 |) \D n+2 ^\ \D n ^\ 

<aj e£+ 2) \D n+2 ^\ 2 J ei n + + 2) (|V^ n+1 , n+2 | 2 + |V 2 ^ n+1 , n+2 |) 2 \D n ^\ 



i rs{n+2) | n ,,2 , -1 2 [ rs{n+2) X{\ x n+1 - X n+2 \ > £) . |2 

< a I e' n+1 \D n+2 iP\ +Ca a 6^ +1 j- — ^ \D n i(>\ 

l^r4 / o( n + 2 ) 



(9.48) 



+ c a - L N" / x (kn+i - x n+2 | < £) |a^|' 



where we used the bounds for |Vw| and |V 2 w| from ()5.4p and that w(x) = a/\x\ for |x| > £ since 
£ S> i?/iV. Using (|9.46p to bound the last term, we obtain 



/ ei n + 2) (|v Wn+1 , n+2 | 2 + |v 2 W n +l,n+2 



1 



(9.49) 



To bound the second term on the r.h.s., we apply Hardy inequality. We have 



p,(n+2) 



1 



\x n+ i-x n+2 \ 2 



Vn + i ( e(«+ 2 ) 



<c| Q^\D n+1 i,\ 2 + C£- 2 f 



\D n ^\ 2 

' 9 n+l ™ \ 2 



8=1 



(9.50) 



< C / e£*+ 2 > \D n+1 4>\ 2 + C(iV - n)- l t 



yn+1 



j>n+l i=l 
2 



<c| e^i^wi 2 + c(N-n)~ 1 r 2 J ei n i|A^| 

Since iW 2 > 1, it follows from (|9~i9"j) that 

©i+ + i 2) (|Vu; n+1 , n+2 | 2 + |V 2 u; n+1 , n+2 |) \D n+2 ^\ \D n ^\ 

<o(i)[J ® { ::?\D n ^\ 2 + f e(r +1) I AmVf + / eSi^l 2 } . 
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Part i) of Lemma 19.41 follows now from (|9.44j) . (|9.47p and from last equation. 
In order to prove part ii) we rewrite the l.h.s. of (|9,42p as follows. 

(1 - w n+ i in+2 ) 2 @! n+2) V N (x n+2 - x n+3 ) \D n+ i(j> n+1>n+2 \ 2 



(9.51) 



> / (1 - w n+ i, n+2 ) & n+ i V N (x n+2 - x n+3 ) \D n+1 (j) n+1>n+2 \ 

Using 



n 1 n / i Vw n+l n+ 2 

l~Wn+l,n+2 {l-Wn+l,n+2) 

we find 



(1 - w n+ i in+2 ) 2 9^ n+2) V N (x n+2 - x n+3 ) \D n+1 cj) n+ltn 



+2 1 



2 



> (1 - a) f ei n + + 2) V N {x n+2 - x n+3 ) \D n+1 iP\ 2 (9.52) 

- Cor 1 J e££ 2) \v Wn+w \ 2 v N ( Xn+ 2 - x n+3 ) \D n ^\ 2 . 

The last term can be controlled by using (|5.4p and that w(x) = a/\x\ for \x\ > I S> R/N by 
J S^^IV^n+i^l 2 V N (x n+2 - x n+3 ) \D n 1p\ 2 

< CN 2 [ e ( ^ 2) x(\x n+1 - x n+2 \ < I) V N (x n+2 - x n+3 ) \D n ^\ 2 



. n 2 f ^(n+2) X{\Xn+l ~ X n+2 \ > I) 2 

+ Ca / 8; +1 : ^ V N {x n+2 - x n+3 ) \D n ip\ 

J |^n+l ^n+2| 

< CN 2 e~ C£ ~ e J V N (x n+2 - x n+3 ) \D r M 2 + Ca 2 r 4 J V N (x n+2 - x n+3 ) \D n ^ 

< o(l) / effi V N (x n+2 - x n+3 ) \D n i>\ 2 . 



(9.53) 

From (|9.5ip . we have 

(1 - w n+ i in+2 ) 2 ©l n+2) V N (x n+2 - x n+3 ) \D n+1 (p n+ltn+2 \ 2 

> (1 - o(l)) / ei n + + 2) V N (x n+2 - x n+3 ) \D n+l ^\ 2 (9.54) 



(1) / e^ +1 ) V N (x n+1 - x n+2 ) \D n i>\ 



2 



In the last term we used B^V^ < @n +1 ^ , the permutation symmetry of ifi and we shifted the indices 
n + 2,n + 3— >n + l,n + 2 □ 

Proposition 9.5. Suppose Ni 2 3> 1. Let {)- n ^ be defined as in \9. 1J$ . Then, if N is large enough 
(depending on n), 



C£N n J @ n +2) D ^i n) 1> D nhf ] ^ + h.e.> C%N n+ \l - o(l)) J Q^ 2) |Vi D n+1 </f 

- fintyO 



i<n<j 



(9.55) 
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where the error term £l n (ip) has been defined in $9.18\) . 
Proof. We rewrite the l.h.s. of (|9.55p as 



i<n<j 

C£N n j ^ { n +2) D n ^i> 

i<n<j 

E J2 J @ n +2) Vn ( X 3 ~ O AA ? D n ^ 



C'$N r > 



2 

i<n<j m>n, m^=j 



- C%N n Y, E A - / ™ n+2) D n(VN(xi - x r ) VO D n Aj ^ 

i<n<j rj^i 

+ YY. X - E J & n +2) A (V N ( Xi - X r ) D n (V N ( Xj - X m ) + h.C 



i<n<j r=£i m>n,m=£j ' 



(9.56) 



with A r = 1 if r > n, and A r = 1/2 if r < n (recall the definition of t)\ n \ for i < n, in (|9.14p ). The 
terms on the last two lines are easy to bound because the potential V/v(xi — x r ) forces the particle i 
to be close (on the length scale iV -1 ) to the particle r. But then the factor #j n+2 ) in Gn™ +2 ^ makes 
this contribution exponentially small. More precisely, for i < n, we have the bound 

(v Q G^ +2) ) |V%(a?i - x r )\ < e- ce ~ s &fr+V (9.57) 

for a = 0,1, (3 = 0,1,2, and for all N large enough. It is therefore easy to prove that 

C%N" Y [ Q n +2) Af^A^ + h.C. 

i<n<j 

c%N n Y j @( n +2) AA^AA^ 

i<n<j 



CJ$N T 



2 

i<n<j m>n, rnj^j 



o (e- c ^) f {e^iA+i^l 2 + e£ilAVf + e^l A-iVf + eir 3 2) |A-^l 2 } 



Lemma 19.51 now follows from Lemma 19.61 and Lemma 19.71 below, where we handle the first and, 
respectively, the second term on the r.h.s. of the last equation. □ 

Lemma 9.6. Suppose the assumptions of Lemma \9.5\ are satisfied. Then we have 

CoN n Y J @ n +2) AA ? AA i^ + h.c. 

i<n<j 

> c %N n+ \i - (i)) J ei n + 2) iVxA+i VI 2 ( 9 - 58 ) 

-o{N n+2 ) [ Q^\D n+1 ^\ 2 - o{N n+l ) / e^iViAVf- 
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Proof. Integration by parts leads to 

E J @ ( n +2) D n Ai i> D n Aj x{; + h.c. 

i<n<j 

= E / e « n+2) |V i V J -£>„i/'| 2 + E [Vie { n +2) ViVjD n ^VjD n i> 

i<n<j i<n<j 

+ E [Vj&n^ViDn^ViVjDnlP 

i<n<j 

+ E f ViVj-e^+ 2 ) ViA^Vj-A^ + kc. 

i<n<j 

The second term on the r.h.s. of the last equation can be bounded by 

i<n<j 

^« E / |v l5 + + 2) )|2 iVjAVf+tt' 1 E /©^IV^AVI 2 

i<n<j ®n i<n<j 

for some q > 0. Next we use that, by Lemma EH part iv), 

|V7.p)( n + 2 )|2 

E 7^ 1 < cr 2 e(r+ 1 ) 

^ o( n + 2 ) 

and therefore, since N£ 2 3> 1, 

E | / v * n" +2) ViVjD^VjD^ 

i<n<j 

<aCr 2 Y [ ® ( n +1) iVPn^ + a" 1 E / ^ +2) IV.V.AVI 2 

j>n i<n<j 

< o(n 2 ) f e^T 1 ) |A+i^| 2 + o(i) e / ^ n+2) Iv^-a^I 2 - 

i<n<j ^ 

The estimate of the third term on the r.h.s. of (|9.59|) is almost identical to the second term 



i<n<j 



^«E /^W^i v ^i 2 +«" 1 E f @i n +2) iv^-a^p 

i<n<j ®n i<n<j 

<c«r 2 E / e^iv.A^P + a' 1 E / - n+2) IV^-AVI 2 

<o(iV) | G^ +1 ) |V!A^| 2 + o(l) E / 6 « n+2) |V,V,AVf • 



i<n<j ' 
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Finally, to bound the fourth term on the r.h.s. of (|9.59p . we use that, by Lemma |A.1| part vi), 

E ' V i Vi " n+2) ! - Cr2 ®n +1) and E l V J' Vl @ n +2) \ ^ Cl- 2 e£ +1 \ (9.63) 

j>n i<n 

This implies that 

Y, | J v J v i e^v l D n lpv j D r ^ 

^E/ El^ V i « n+2) H V ^^| 2 + E / El V ' V ^ « n+2) )l l V i^Vf (9.64) 

i<n j>n j>n i<n 

<o(N) [ e( n+1 )|ViA^| 2 + o(iV 2 ) / G^lDn+i^l 2 . 



i<n<j 



Lemma EH now follows from (19391) . (|9Uil . (19U2D and (IQ4"]) . □ 

Lemma 9.7. Suppose the assumptions of Lemma \9.5\ are satisfied. Then we have, for N large enough 
( depending on n ), 



C%N n 



E E J @ n +2) Vn{ - X 3 ~ X ^ D ^ ^D n lP + h.C. 
i<n<j m>n, mj^j 

> -o(N n+3 ) J ®^V N (x n+1 -x n+2 )\D n ^\ 2 . 



i<n<j m>n, mjtj ' 



Proof. We have 

E E J Q n +2) V n(Xj ~ *m) Dn&i ? D n $ + h.C. 

= E E / 9 ^ +2) V N (x m - Xj )\ViD n ip\ 2 (966) 

i<n<j m>n, m^j 

+ E E / V * " n+2) VN ^ ~ Xm ) ViDn ^ D n^ + h-C • 



i<n<j m>n, m=£j 

The second term can be bounded by 



i<n<j m>n, m^j 



E E / v * - n+2) Vn ^ - x ™) VA ^ D » ^ + hx - 

<« E E [ lV ^ 12 V N ( Xj - Xm )\D n ^ (9.67) 



i<n<j m>n, m=£j ' 



i<n<j m>n, m^j 

Since, by Lemma lA.ll part iv), 



+ a ^ E E / - n+2) v *( x i - x -) i V * D » 
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using the permutation symmetry and optimizing a, we obtain 

I E E / Vi 9 « n+2) " Xm) ^* D n^D n i> + h.C. 

< o(N 3 ) / e(r +1) V N (x n+1 - x n+2 ) \D n ip\ 2 



i<n<j m>n,m^j ' 



(9.69) 



+ E E /e^^^-x^iv^vi 2 - 



i<n<j m>n, m^j ' 

Inserting the last bound in (19.66f) . we conclude the proof of Lemma 19.71 



□ 



A Properties of the cutoff function 9 



(n) 



Recall that the cutoff functions 6^ n) = 6^ n) (x) defined for k = 1, . . . , N and n € N, in Eq. (|532l) . 
In the following lemma we collect some of their important properties which were used in the energy 
estimate, Proposition 15.31 



Lemma A.l. i) The functions @^ are monotonic in both parameters, that is for any n, k G N ; 



©Si < ©i n) < i 



e 



(n+l) 



< G<f } < 1 . 



(n) 

Moreover, Q k is permutation symmetric in the first k and the last N — k variables, 
ii) We have, for any n G N, k = 1, . . . , N, 

k N 



EE".. 



e<"> < c m e'"-" 



i=l j^i 

Hi) For every k = 1, . . . , N, n S N, we have 



\^& { k n) \<ct-U^h n )e^<ct-^ 



2" 



N 



-1q("-1) 



r=l 



(on ^ \ 
tfE^J 9 



if i < k 
if i > k 



iv) For every k = 1, . . . , N , n E N we /iav e 



iv | Vi G fc 



(n) 



E 



For every fixed k = 1, . . . , N and n £ N we /ictue 



v 4 v,ei n) 



-2 



5> 



2 n 



m=l 



/iri ©i n) < CI' 2 6^ n 1} , if i^j and i,j > k 



(A.l) 



(A.2) 



(A.3) 



r=l 



< cr z e 



2 O^" 1 ) 



for any i,j 



(A.4) 
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vi) For every fixed k = 1, . . . , N and neNrae have 

E h v ; 



(A.5) 



Proof. Part i) follows trivially from the definition of 0^ n \ Part ii) follows from x m e x < C m e x l 2 
for every real x. To prove part iii), we observe that, for i > k 



h Ik 

2 n \ / 2 n 



r=l j^r 



r=l 



Since |V/i(a;)| < 1 h(x), we obtain 



< 



Similarly, for i < k, we have 



v-e (n) 



(n« ^ \ / on ™ 

f Em ex p - 
r=l / \ r=l j^r 



— E + r?r) exp — - EE h 



r=l 



r=l jr'^r 



with ry r = if r > and rj r = 1 if r < A;. Therefore, in this case 

iV \ / „„ AT 



< 



exp 



r=l 



E E v 

r=l j^r 



Eqs. (|A.7j) and (jA.9l) . together with part ii), prove (IA.2j) . 
As for part iv), we have, from (|A.7p . 



N V -P)( n ) 
— l v J u fc 



E 



A? 



j=fc+l K ~> k 



< 



c- 2 E fEv 

j=k+l \ r=l 
( N k 

- or * I F S Ev 

< cr 2 e< n_1) 



) exp | 




I exp 


/ nn ™ \ 



(A.6) 



(A.7) 



(A. 



(A.9) 



(A.10) 



by part ii) of this lemma. The contribution to (|A.3p from terms with j < k can be controlled 
similarly, using (|A.9P • The proof of part v) and vi) is based on simple explicit computations and the 
same bounds used for part iii) and iv). □ 



B Example of an Initial Data 

In this section, we denote by (1 — uj(x)) the ground state solution of the Neumann problem 



-A + -V N {x) ) (1 - co(x)) = e/ (l - u{x)) 
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on the ball {x : \x\ < £} with the normalization condition uo(x) = if \x\ = I. We extend lv{x) = 
for all x G M 3 with |x| > £. We will choose £ such that a <C £ <C 1. Recall that a = a$/N is the 
scattering length of the potential V^(x) = N 2 V(Nx). Assuming that V > is smooth spherical 
symmetric and compactly supported, we have, from Lemma A. 2 in [8], the following properties of ee 
and oj(x). 



i) If a/£ is small enough, then 

ii) There exists cq > such that 
for all i6K 3 . Moreover 

< Ca 



ei = 3ar 3 (l + o(a/£)) 
Co < 1 — uj{x) < 1 
l(\x\ < I) 



\x\ + a 



and \Vuj(x)\ < Ca 



l(\x\ <£) 
(\x\ + a) 2 



N 



i<j 



We define the TV-body wave function 

Wat(x) := - u(xi - ay)) . 

For m = 1, . . . , N, we also define 

W^\x m+ x, . . . ,x N ) := ] [ (l-u(xi-Xj)). 

Lemma B.l. Define 



N 



m<i<j 



\\w N ^)U"LM x 



^at(x) 



for any <p G -fT 1 (M 3 ) with \\<p\\l2 = 1. Then, i/a <C f < 1, we /iat>e 

{^ N ,H N ip N ) <CN 

and, for any fixed k, 



lim U N -<p® k ®^- k) \ 

N— >oo 



o, 



where 



(JV-fc) 
AT ' 



Sfc+i, . . . ,xatJ := 



rifc<^ <J ( i - ^Qg« - gj)) n^fc+i <p{ 
i n*<i<i(i - ^ - x j)) nf= fc +i 



(B.l) 



(B.2) 



(B.3) 
(B.4) 



Proof. Let 0at(x) := Il^i ^(^i). and > for m = 1, . . . , iV, 0^ (x m +i, . . . , x N ) := Hf >m tp(xj)- We 
start by noticing that 



(l-o(l)) 



W 



[i] M 

N 



< 



W N < 



< 



w. 



(B.5) 
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Here IIW^Vjv H ^ s the norm on L 2 



The upper bound in (|B.5P is clear since 1 — u < 1 and 



|c^|| = 1. To prove the lower bound, we note that, by (|B.2j) . and using the notation lu,, — ^i.r, — .r 

TV 

||Wjvcrf = / dx J](l-^-) 2 |^(x)| 2 

i<j 

TV / TV \ TV 

= / dx n (i-^i^x)! 2 - /dx i-na-^j) a IK 1 

Ki<j J \ j=2 / Ki<j 



TV 



3=1 



> M\ 2 \\W l ^^\\ 2 -2J2 / dx Wlj - ^ 1J (x 2 ,...,x w ) |^(x) 



> II^JvVSv 1 !! 2 - CNa / dx 



l(|xi - xj\ < £) 



W$\x 2 ,...,x N j\ |</>tv(x)| 2 



WJ1]||2 



>{l-CNatMUW y ^ 



using that l(\x\ — Xj\ < t) < l\x\ — xA 1 , and then applying a Hardy inequality in the variable x\. 
This proves (|B.5|) . because I <C 1. Analogously, we can prove that 



[fc] Ak] 



TV 



< 



w N 4> N 



< 



w [k] M 
w N q> N 



(B.6) 



where Ofc(l) — > as tV ^ oo, for every fixed k > 1, and where ||W^^S || is the norm on L 2 (R 3 ^ 
Next we prove (|B.4p . To this end we remark that, by (|B.6|) . 



W N 4> N W n 4>n 



\\Wn4>n\ 

as N — > oo. Moreover, since 



\W 



[k] 



TV 



< 



\\W N cP 



TV 



\w. 



[k] 



TV 



(B.7) 



\w 



TV 



we observe from ()B.7p and (|B.6[) that 



limsup ||^)jv — tp® k <8> ^ || 2 < limsup 

TV^oo TV-+00 



Now we have 



(i^v - w; 



Jf ] )0iv 



dx II- [] (1 

i<j<k,i<k<j 

N 



\W [k] 6 [k] \\ 2 
\ VV N II 



[^if ] (x fe+1 ,...,x JV )] 2 ni¥'(^)i s 



(Bi 



v 



13) 



/V . TV 

<c^2Y^ / dx^fwf^xfc+i,...,^)] 2 !!^^)! 

i<fe i=l J j=l 



i<k j=l 

<CNka£y\\ 2 Hl \\wP^ ] \\ 2 



by using (|B.2p and Sobolev inequality in x\ (see Lemma f6.4l part i)). By (1B.8P and I « 1, this 
proves (|B.4p . 



1 



Finally, we prove (|B.3p . To this end we observe that 

N N 



N N 



w N 

where 



h n (w n 4> n ) = ^2L j 4> N + e e '^2i(\x m -x j \ <e)4> N -^2 E 7 

1=1 j^m^i^^m J 



m£3 



Note that 



w^4> N Ljif; N = I wjj Lj4> N ii) N = I w^Vj4> N VjipN- 



From (|B.9p we find, by using (jB.ip . W^r < Wjy and by applying the Sobolev type inequalities of 
Lemma 16.41 and the permutational symmetries, 

(W N 4> N ,H N W N 4> N ) 

N . N . 

tv| 2 + ee £ / dx^(x)l(|x j -x m \ < 



N 



(x) 



N N 

E E 

»=1 j,"i^i,j^m 



i — • 1 

1 — Wjj 1 — Wjrn 



+ CiV(iV- l)a|| (p\\%i 



+ CN(N - 1)(N - 2)a 2 \\ipf Hl 
for any e > 0. From (|B.6p . and since f < 1, we have 



<V 3] 



||Wjv0jv| 



IIW^jvII 



which completes the proof of (|B.3[) . 



(B.10) 



(B.ll) 
□ 



C Trapped condensates 

In this Appendix we show that Theorem 12.21 can be applied to the ground state of interacting Bose 
Hamiltonians with a trap. Recall the definition of the Hamiltonian Hn without a trap from (|2.ip . 
and define 



N 



N 



N 



j=l j=l i<j 

with a smooth trapping potential V cx t > satisfying limi a .i_ >00 V ex t(x) = oo . Denote by the pos- 
itive normalized ground state vector of ff^ ap . The corresponding Gross-Pitaevskii energy functional 
is given by 

4p P (£>= / dx (|V < /.(x)| 2 + yext(x)|0(x)| 2 +47rao| < /.(x)| 4 ) 
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and we denote by (frcp the L 2 -normalized, positive minimizer of £ G r p P . As proven in |16| . the ground 
state energy per particle is given by minimum value of <?Qp P as N — ► oo, 

l(^ aP , H^T) - 4TOgp P )> (C.12) 

and the one-particle marginal density 7jJ trap associated with Vv aP satisfies 7^ trap — > |<^gp P )(0gp P | 
(with convergence in the trace-norm). From (|C.12|) . (Vv aP > H^^tjj^) < CN and since Hn < H^ av , 
we obtain that Vv ap satisfies (|2.17p . The goal of this section is to prove in Proposition IC.2I below 
that V^ ap satisfies the asymptotic factorization property (|2.18j> . From Theorem 12.21 we therefore 
immediately obtain the following corollary: 

Corollary C.l. Suppose V satisfies the same conditions as in Theorem \2.S[ LetipN,t be the solution 
of the Schrodinger equation without a trap, idtipN,t = HNipN,t> but with initial data given by the 
trapped ground state, ipN,o '■= ip% ap ■ For k = 1, . . . ,N, let 7^ be the one-particle marginal density 
associated with tpN,t- Then, for every t 6l, and k > 1, 

7$-|Pt>W fc as 00 (C.13) 
in the weak* topology of C 1 (L 2 (M. 3k )) . Here (ft is the solution to the Gross- Pitaevskii equation 

id t <pt = -A<ft + 8na \(pt\ 2 (pt 

with initial data (ft=o = ^gp ■ d 

Proposition C.2. For any fixed k = 1, 2, . . there exists a sequence of normalized wave functions, 
^~ k) G L 2 (R 3( - N ~ k )), N > k, such that 



\*T-[fig\* h ®<fN~ k) \ 



as N — > 00. 



We will prove this proposition only for k = 1, the proof for arbitrary k > 1 can be obtained 
similarly. For brevity, we set £jy = £/v _1 - ^ or the proof, we make use of the following three lemmas. 

Lemma C.3. There exists a constant C > independent of R,N such that 

\\l(\ Xl \ > R)^% ap \\ < Ce- R (C.14) 

where l(s > A) denotes the characteristic function of the interval [A, 00). 

Lemma C.4. We have (x) > for all x £ M 3 . Moreover 

||(1 - A)0*™f|| < 00, i^cp, V ext {x)^%f) < 00 
and there exists a constant C > such that 

\\l(\x\ > < Ce~ R 

for all R>0. 
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Lemma C.5. For fixed R > 0, N e N define £ R<N G L^R 3 ^" 1 )) 6y 

1 



£,r,n(x-n- 



where xjv-i = (^2> • • • ,xn). Then we have 



/ da?i |0 G /(xi)| 

7|xi| 



dx 



iV-1 



cLei 



|xi|<i? 



^t? p (a;i,xjv-i) - ^Sp { x i)£r,n(*n-i) 



< c R d N 



(C.15) 



where c R < oo is independent of N and d^ is independent of R and satisfies d^ — > as N — > oo 

Using these three lemmas we can prove Proposition IC.21 
Proof of Proposition \C.S\ for k = 1 . Using the notation introduced in Lemma IC5I we have 

dxi |V>^ aP (>i,Xjv_i) - (/)Qp P (a;i)^ Ri Ar(xAr_i)| 2 

dxi |Vw aP (x) - </>%'p(xi)Zr ) n(x.N-i)\ 2 
, trap^.^ ^ ra P 



dxAT_i 
dXjV-1 
- f dxAT_i 

<c R d N + Ce~ R 



\xi\<R 



/ dxjv_i / dxi |^ ap (x) - ^Gp P (zi)^jv(xAr-i)|' 

J J\ Xl \>R 



(C.16) 



where we used Lemma [C.5l to bound the term on the second line, and Lemmas IC. 31 and IC.4I to bound 
the term on the third line. Eq. (|C.16p implies that 



£r,n 



R,N\ 



,/rap 



< 



^ap 



i>Gp P ® £fl,iv|| 
0GP P ® £r,JV 



(C.17) 



1 



^trap 



Now choose a sequence .Rat such that i?7v — > oo and c RN dN -> as iV -> oo. Then, taking 
£iv = C.Rjv,iv/||^iijv,iv||; we clearly have ||£jv|| = 1 for all N, and, by (|C.16jl and <[C.17|) . 

</w p <2D ^"iv 1 1 — ^ as N — > oo . 



& trap 



□ 



We still have to prove Lemmas IC.31 IC.4I and IC.5L Lemma IC.4I is a standard result which follows 
from the fact that </>Qp p is the solution of the elliptic non-linear eigenvalue equation 



a /trap . t t /trap . I j^rapiz /rrap /rrai 

A< Pgp + Kixt^QP + 8vra o |0 G p p | GP P = ^0 GP ' 



.trap 



trap 1 2 jtrap 



,trap 



(C.18) 



with some constant [i. Lemma IC.5I has been proven in [16], more precisely, it follows from Eq. (13) 
of [16] by noticing that the two terms in the parenthesis in this equation converge to zero, uniformly 
in R, because of Eq. (7) and Lemma 1 in [16J. It only remains to prove Lemma IC. 31 To this end we 
use the following two lemmas. 
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Lemma C.6. Let x G C°°(R) wita x(s) = if s < 1 and x(s) = I if s > 2, and let f E C^R) 6e a 
monotonically increasing function with sup x \ f'(x)\ < oo. T/ien we Ziave, /or R > Zarye enough, 

X (\xx\/R) (< rap - |/'(M)| 2 - x(|*i|/22) > x(N/i?) 2 , 



where En denotes the ground state energy of Ht Tap 



Proof. Define 

N N 

= (~ A i + Vextfa)) + V N (Xi- Xj ) 
j=2 2<i<j 

and let En-i = inf a(H^^,). Moreover, we define ipj^i & L 2 (M^ N ~^) to be the positive normalized 
ground state of H^^. Then we have, since — Ai > and Vn(x) > 0, 



Xdx^/R^H^-If'dx^-E^x^l/R) 

> X (\xi\/R) (H^ + V ext ( Xl ) - \f(\ Xl \)\ 2 - E N ) X (\xi\/R) (C.19) 

> x(\xi\/R) 2 (Vcxt(zi) ~C-{E N - E N -i 



where we used the assumption \f'\ < C. Next we remark that there exists a constant C > such 
that 

E N < E N -i + C for all N . 
In fact (using the symmetry of the wave function) 

E N < (<fc? ® Pn%Hn <t>%7 ® fc) = En-i + ("Ai + V c Axi)) *gP) 

+ (C p ® (JV - 1)JVV(% - x 2 ))^ ® fe) 

< + C ||(1 - A)^ p || 2 + C(<^f , y e xt(x!)^p P ) 

< ^-i + C 

(C.20) 

where we used the operator inequality W{x\ — x%) < C||W||ii(l — Ai) 2 and Lemma [C.4I Since 
limij.i^jx, V ex t(x) = oo, the lemma now follows from ()C19|) . □ 

Lemma C.7. Suppose that f,x o,re as in Lemma \C.(A Then we have, for R large enough, 

\\e^x(\xi\/R)^ ap \\ < C R (C.21) 

for some constant Cr depending on R but not on N . 

Proof. We compute 

^(M^trap _ EN)e -f(\ Xl \) = tftrap _ | /( | Xl |)|2 — En + i (pi ' ^r/'CH) + AM A • Pi) , 

V l^i I l^il / 

with pi = — iSJ\. Therefore, for R large enough, 

Re(ef^ X (\xi\/R)^T^ filXll) ( H T P ~ e n) e~f^ef^x(\ Xl \/R)^) 

= (e^H% ap x(\xi\/R), (h^ p ~ \f'(\xi\)\ 2 - E N ) X (\x 1 \/R)ef^^) 
> \\ef^x(\xi\/R)i> N ap f 

(C.22) 
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where we used Lemma IC.6I On the other hand 

Re(e^ X (\ Xl \/R)^ p , [h^ p - E N ) e~^ e^x{\ Xl \/ R)^ p ) 



< \\e f ^x{\xi\/R)ip 



trap I 

AT 



H%^, X {\ X1 \/R) 



because (H^ p - E N )^ P = 0. Combining (1C321) and (RT25I) we obtain that, for R large enough, 



^M) x (\ Xl \/R)ifi 



trap I 
N 



< 



H N , X (\xi\/R) 



trap 
N 



Next we note that 



[H N ,x{\xx\/R)\ = -2iR~ x x!(\xi\/R)p-\ ■ V l +R~ 2 x"(\x 1 \/R)+R- 1 ^^- 

\xi\ \x\\ 



Since / is monotone increasing, we see that 



^ X '(\xi\/R) 



\ x i\ 



<Ce^ 2R \ 



J(\ X1 \) X'(\xi\/R) 
\ x i\ 



< CR~ l e f{2R) and 



Jilxil) x"(\xi\/R) < Ce f(2R l 



(C.24) 



The energy estimate (|C.12|) and Vn > imply that ||Vi?/>jv aP |l < C uniformly in N. From these 



estimates the lemma follows. 



□ 



Proof of Lemma \C.3l Suppose x is as m Lemmas IC.6I and IC.7I For a fixed Rq large enough, we have, 
by Lemma IC.7} 



le^l^H < \\e^x(\xi\/Ro)^ p \\ + He^l (1 - x(\xi\/R )) 4>7 P \\ < & 



Therefore 



l(|xi| > «)Vw aP H < ||e" |a;i| l(|xi| > R)e^iP% ap \\ < Ce 



-R 



□ 



D Properties of the one-body scattering solution 1 — w{x) 

In this section we prove part i) and iii) of Lemma 15,11 



Lemma D.l. Suppose that V > is smooth, spherical symmetric with compact support and with 
scattering length do- Let 

poo 

p = suprV(r) + / drrV(r), (D.l) 
r>o Jo 

and suppose (fo(x) is the solution of 



-A + —V ) ifo = with (fo — > 1 as \x\ —> oo . 



(D.2) 
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i) There exists Co > 0, depending on V , such that Cq < <po(x) < 1 for all i£R 3 . Moreover there 
exists a universal constant c such that 



1 - cp < (fo(x) < 1 for all x eR 3 . 
ii) There exists a universal constant c > such that 

|V^o(^)| < C-j— \V ipo(x)\ < C-^-r and \V 2 ipo(x)\ < C-r^pr 

\x\ z \x\ \x\ A 

Moreover there are constant C\,Ci, depending on the potential V, such that 

\V<po(x)\ < Ci |vVo| < C 2 . 



(D.3) 



(D.4) 



(D.5) 



Proof. Let R be such that suppF C {x 6 M 3 : \x\ < R}, and let ao denote the scattering length of 
V. Then we fix R > R such that glq/R < min (p, 1/2), with p defined in (fDTTT) . 

In order to prove part i), we observe that, for |x| > R, (po( x ) = 1 — oo/l^l- Hence 

1 



< (fo(x) < 1, 



and 



1 - P < fo(x) < 1 



for Ixl > R. 



(D.6) 



Next, by Harnack principle the ratio between the supremum and the infimum of <po in a given ball 
is bounded: therefore <po is bounded away from zero in the ball \x\ < R and thus there exists 
Cq > such that 9?o( x ) > Co f° r an x £ Moreover by the maximum principle, and since, from 
pX2]h -A<p < 0, it follows that <p (x) < 1, for all x £ M 3 . To prove ([DT3]) for |x| < we write 
fo(x) = m(r)/r, with r = |z|. Then m'(R) = 1, and, from (jD.2j) . 



— m" (r) + — V(r)m(r) = . 
Since < <po(x) < 1, it follows that m(0) = and < m{r)/r < 1. Therefore, for r < R, 



(D.7) 



m'(r) = m'(R) — / dsm"(s) = 1 



dss y( s )!^£2 > i-c 



ds s V(s) > 1 - cp (Di 



and 



m(r) 



dsm'(s) > r(l — cp) 



?"(r) > 1 



cp 



for all r < i? . 



(D.9) 



The last equation, together with ()D.6|) . implies (|D.3|) . 

Next we prove ii). For \x\ > R, we have <po(x) = 1 — ao/M and thus 

|V<po(^)| < < < A, for \x\ > R, 
fI R\x\ \x\ 

by definition of R. Next, for \x\ < R, we write <po( x ) = rn(r)/r, with r = |x|. Then 

m'(r)r — m(r) 



(D.10) 



|V<po(*)| 



— asm is) =■ / as aum (k) 

r Jo r Jo Jo 

1 r 

~2 / d/t/sm (k) , 
r Jo 



(D.ll) 
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because m" '(k) > 0. From (|D.7p we obtain 



|V^ (x)|<^^ r d 



a 



KK V(k) < C- — 



(D.12) 



because 87rao = f V(x)ipo(x) (see Lemma l5.ip . part iv). Moreover, again from (jD.lip and (|D.7p . we 
have 



\V<Po(x)\ < 



1 



2r 2 Jo 



dK K 2 V(k) 



m(«) ^ sup k > k 2 1/(k) p 



< c 



< c- 



(D.13) 



Together with (ID.lOh we obtain the first two inequalities in (ID,4j) . From (|D.10p and from the first 
inequality in (|P.13p . it also follows that there exists C\, depending on the bounded potential V, such 
that |Vyo(aO| < Ci- To prove the second bounds in (ID.4|) and (1D.5|) . we note that 



IV72 f \\ ^ a o . P 

|v v?o(x)l -w w 



for Ixl > i?, 



(D.14) 



by the definition of R. For \x\ < R, we have (expanding m{r) and m!(r) and using that m(0) = 0) 

|VVo(s)| < 



m"(r) m'(r) m(r) 
2 ^ "t - 2 ^ 



2 w r 



2 r 

+ ^ / dss 2 F(s 
^ 3 Jo 



m{s) 



(D.15) 



, (sup s>0 s 2 ^(s)) ^ /9 

< C =— 7t < C-7T 



Last equation, together with ()D.14p . implies the third bound in (|D.4p , Moreover, from ()D.14p and 
the second line in (ID.15j) . it also follows that there exists C2, depending on the bounded potential 
V, such that \V 2 (p (x)\ < C 2 . □ 

Proof of Lemma \5.1\ part i) and in). By scaling l — w(x) = (fo(Nx), with ipo defined in Lemma [D. II 
Therefore part i) of Lemma [5.11 follows immediately by part i) of Lemma [D. 11 and part iii) of Lemma 
ED follows from pTIj) and (|L\5|) . □ 
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